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ABSTRACT 



Aims. We present new VUT ISAAC spectra for 30 quasars, which we combine with previous data to yield a sample of 53 intermediate 
redshift (z ~ 0.9 - 3.0) sources. The sample is used to explore properties of prominent lines in the Hj8 spectral region of these very 
luminous quasars. 

Methods. We compare this data with two large low redshift (z <0.8) samples in a search for trends over almost 6dex in source 
luminosity. 

Results. We find two major trends: (1) a systematic increase of minimum FWHM H/? with luminosity (discussed in a previous paper). 
This lower FWHM envelope is best fit assuming that the narrowest sources radiate near the Eddington limit, show line emission 
from a virialized cloud distribution, and obey a well defined broad line region size vs. luminosity relation. (2) A systematic decrease 
of equivalent width [Oin]/l/i4959,5007 (from W»15 to ~1A) with increasing source bolometric luminosity (from logLboi ~ 43 to 
logZ-toi ~ 49). Further identified trends require discrimination between so-called Population A and B sources. We generate median 
composite spectra in six luminosity bins to maximize S/N. Pop. A sources show reasonably symmetric Lorentzian HyS profiles at all 
luminosities while Pop. B sources require two component fits involving an unshifted broad and a redshifted very broad component. 
Very broad IIy6 increases in strength with increasing log L|,oi while the broad component remains constant resulting in an apparent 
"Baldwin effect" with equivalent width decreasing from W ~80 to ~20 A over our sample luminosity range. The roughly constant 
equivalent width shown by the IIy6 very broad component implies production in optically-thick, photoionized gas. The onset of the 
redshifted very broad component appears to be a critical change that occurs near the Pop. A-B boundary at FWHM Hfi 4000 
km s"' which we relate to a critical Eddington ratio (« 0.2±0. 1). 

Key words, quasars: general quasars: emission lines - emission lines: profiles - black hole physics 



1. Introduction 

Much remains to be defined in quasar research even if one re- 
stricts attention to the broad emission hnes often used to de- 
fine them. Their broad line spectra show considerable diversity 
complicating attempts to generate composite spectra and mak- 
ing estimates of intrinsic properties such as black hole mass un- 
reliable. In addition, the relationship between broad line emit- 
ting active galactic nuclei (Type 1 AGN) and various classes of 
sources that do not show broad lines (e.g. Type 2 AGN, LINERs, 
Blazars, NLRGs) is still uncertain. It is difficult to imagine how 
advances in physical understanding can come until source phe- 
nomenology is clarified. We have focused on clarifying the phe- 
nomenology of Type 1 sources because of their relatively unam- 
biguous broad line signature. Leaving aside reverberation stud- 
ies making use of source variability, observational advances can 
be expected to come from two areas of spectroscopic investi- 



gation: (1) inter-line comparisons (from Lyg to Hg) i i i terni s 
of relative intensity and profile shape (e.g., IShang et"an , l2007h : 
(2) moderate dispersion and high S/N studies of the Hfi spec- 
tral region in quasars with z > 1.0. There is a long history 
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cently has it become possible to obtain spectra with resolution 
and S/N comparable to those for optical spectra of low z sources. 
IR spectroscopy of H/3 at z ^1.0 enables one to use the same line 
(H^), the same rest frame determination ([Oiii]/l/l4959,5007 or 
narrow component of li/3) and the same reduction procedures for 
estimating the central black hole mass (Mbh)- Even so observa- 
tions of the most numerous and most luminous quasars at z ~ 2 
can only be made when the optical lines are redshifted into one 
of the IR bands of high atmospheric transmission. 

Spectrographs on HST have provided UV coverage of higher 
ionization broad lines in low-z Type 1 sources (Seyfert 1 
and quasars) with ~150 sources having useful measures of 
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CiV/11549 (e.g iBrotherton et al.L [19941 iMarziani et al.L 119961: 
iBaskin & Laoii 120051: ISulentic et aLl l2007h . Almost all of these 
sources now have mat ching optical spectroscopic measures of 
the HB region (e . g.. [ Marziani et aU, 12003 at iBaskin & Laoit 
l2005l:IShang et al.Ll2007l) . Continuing with a Htany of difficulties 
we note that even some of these UV observations show marginal 
S/N while many do not cover the full range from Lya to 3500 A. 
Optical coverage is usually not synoptical. This makes detailed 
intercomparisons of the strongest high and low ionization broad 
lines possible for only a few tens of objects at best. The lack of 
spectra covering the rest frame from Lya to Ha is unfortunate 
since line intensity and profile ratios provide a wealth of physi- 
cal constraints. The ability to compare lines from ions of widely 
different ionization potentials has helped to elucidate two main 
emitting regions within the broad line region (BLR): (a) one re- 
sponsible for the production of low-ionization lines (LIL) like 
Fen, Hy6 as well as Mgii/12800 and (b) one that emits mainly 
lines (Civ/11549, Heii/t4686) from ions of hi gh ionization po- 
tential (HIL; e.g., IColhn-Souffrin et all [19881) . Meaningful line 
profile studies begin with observations of Civ/I1549and H/3 as 
the most typical HIL and LIL. 

Our attempt at clarifying type 1 AGN phenomenology in- 
volves the Eigenvector 1 formalism o fl oroson & Greed (Il992h 
and was later expanded into 4 dimensions (4DE1). The lat- 
ter f ocused on four param et ers inc luding measures of H/3 and 
Fen ( ISulentic et al.L l2000a|[bl |2007[) . Changes in the width and 
relative strength of li/3 and Fen lines appear to be primar- 
ily r elated to Eddington ratio convolved with source orienta- 



tion jShang et al.Ll2003l: [Marziani et al.L [2003bl:l Yip et"ail 120041: 



iBorosonL \200± icollin et all l2006ir 4DEl studies also intro- 
duced the concept of two quasar populations A and B that 
maximize phenomenological differences and possibly identify 
sources with higher and lower Eddington ratios respectively. 
The width of broad Hy6 (FWHM or 2"'' profile moment cr) have 
been widely used as measures of velocity dispersion in the line 
emitting gas (see e.g. Peterson et al.Ll2004t ISulentic et"an , l2006l : 
IVestergaard & PetersonL l2006t) . allowing what are thought to be 
the most reliable estimates of black hole mass in low z quasars. 
This approach has motivated us to seek similar measures in 
sources with the highest possible redshift rather than using other 
lines as Hfi surrogates. 

In this paper we leave aside inter-line comparisons and fo- 
cus on the H/3 spectral range in low and in intermediate-to- 
high z quasars (z ^ 1). Near IR spectra of the latter were ob- 
tained with the infrared spectrometer ISAAC on ESO VLT Unit 
1. The high S/N and resolution of the new observations al- 
low a meaningful comparison with measures obtained for low 
z sources. Low z measures come from two samp l es: (1) an at- 
las of bright sources with z£ 0.8 (IMarziani et al.L [2003 al here- 
after, ATLAS sample) and (2) a magnitud e-limited sample of 
321 SDSS quasars (mg <17.0, z <0.7 from IZamfiret al.Ll2008l 
hereafter SDSS sample). 

In two previous papers we focused on a search for lu- 
minosity effec t s invo lving H/3, FeUopt, and [Oin]/l/l4959,5007 
(ISulentic et al.L |2004 hereafter Paper 1) a s well as on the use 
of H;8 as a virial black hole mass estimator (ISulentic et al.L[2006[ 
hereafter Paper II). This paper presents observations and reduc- 
tions for 30 additional sources (z - 1.08 - 3.09; We briefly 
discuss the H/3, Fen and narrow line measures (® and then show 
that the addition of 27 new sources to the VLT sample reen- 
forces the luminosity trends described in Paper I (® . We then 
consider the distribution of VLT sources in the optical plane of 
4DE1. Median 4DE1 and luminosity binned composite spectra 
are discussed in the context of the 4DE1 Population A-B concept 



(ISulentic et al.Ll2000al[2008l) . In g5]we use the binned composite 
spectra to decompose H/3 and make estimates of black hole mass 
Mbh and Eddington ratio Lboi/^Edd for Pop. A and B sources. 

2. Observations and Data Reduction 

Spectra for intermediate-high redshift sources were obtained 
between 07/2004 and 07/2006 in service mode with the in- 
frared spectrometer ISAAC mounted on VLTl (ANTU) at the 
European Southern Observatory. Table \T\ summarizes the new 
observations. The basic format is given below the Table, fol- 
lowing almost exactly Tab. 1 of Paper II. In addition to the for- 
mat provided below Table [1] we note that Col. (3) lists the blue 
apparent magnitude t aken from Hamburg- ESO survey papers 
dReimers et all 119961: [wTsotzki et al.L |200(]|) while Col. (4) Hsts 
the source redshift z computed as described in section 12.21 Col. 
(8) lists the radio-loudness parameter /?k defined as the ratio be- 
tween the specific flux at 6 cm and 4400 A in the rest frame. We 
applied a A:-correction for both the radio and optical data. In the 
case of the radio data: /e,6cm = /v.obs ■ [vobs(l + z)/v6cm]" where 
the radio spectral-index is a = 0.5. Only 14 of the sources have 
radio detections; 1.4GHz upper limits were derived for unde- 
tected sources from the NVSS detection thre shold or the Sydne y 
University Molonglo Sky Survey (SUMSS; [Bock et al.L [l999l) . 
Two sources are fo rmally radio-loud following our definitions 
dZamfiretalL 120081) . 

Each spectrum corresponds to a wavelength range (IR win- 
dows sZ, J, sH, sK; Col. (10) of Tab. [T]i that covers all or part 
the region involving HyS Fen/14570 and/or Feii/15130. Reduction 
of quasar spectra and standard stars followed exactly the same 
procedures described in Papers I and II. Wavelength calibration 
yielded rms residuals of 0.4, 0.6 and 0.9 A in the sZ, J and 
sH windows, respectively. Absolute flux measures will be inac- 
curate because atmospheric seeing almost always exceeded the 
slit width (x 0".6) resulting in significant light loss. Small off- 
sets were present in the wavelength calibration because the arc 
lamp frames were obtained in daytime and therefore usually af- 
ter grism movement. A correction for these shifts was obtained 
by measuring the centroids of 2-3 OH sky lines against the arc 
cahbration and calculating the average difference. 

2.1. Sample Considerations 

All 53 sources in the VLT-ISAAC sample were selected from 
the Hamburg-ESO (HE) quasar survey which was a color- 
selected and magnitu de-limited (tob ~ 17.5) quasar survey 
(IWisotzki et al.L [20001) . Thirty-one new spectra are presented for 
30 sources in ^with 27 sources not previously observed. They 
are merged with Paper I and II measures to yield a total sample 
of 53 sources. Two of the new observations involve previously 
observed sources that were discussed in Paper I. The indepen- 
dent pairs of measures provide a valuable consistency check and 
are treated as independent data points in our statistical analysis 
of ^ Two new spectra obtained for HE 1505H-0202 are listed in 
Table[T]with average values used in the analysis. HE 0353-3919 
(Paper I) has been excluded from quantitative analysis because 
its H/? profile is compromised by a gap in coverage between the 
sZ and Z bands (see Fig. 2 of Paper I). The resultant sample of 
52 sources (and 54 data points) is henceforth referred to as the 

ISAAC sample. 

The samples of IMarziani et al] ( 12003 aL ATLAS) and 
IZamfiret all (120081 SDSS) are used as low-z spectroscopic com- 
parison samples. SDSS is more complete, involving the 321 
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brightest g-band selected (<17.0) in SDSS DR5 with z <0.7 
while the ATLAS sample involves a more heterogeneous selec- 
tion of 215 quasars largely brighter than - 16.5. We make 
wide use of the SDSS sample in this study because of its higher 
level of completeness, larger number of sources, and because re- 
sults from the ATLAS sample were already presented in Paper 
I and II. Both the SDSS and ATLAS data enable: (1) accurate 
measurements of Hy6 and Feiiopt emission, (2) decomposition of 
the broad Hfi profile, and (3) measures of the narrow H/? com- 
ponent (HySNc) and [Oiii]/l/l4959,5007 Unes. The SDSS sample 
suff'ers a strong Malmquist bias. The ATLAS sample shows a 
more uniform luminosity distribution making a correlation anal- 
ysis less biased. It will be used in the analysis of luminosity 
trends involving FWHM HyS (©. 

2.2. Data Analysis 

Analysis of ISAAC spectra made use of IRAF tasks to accom- 
plish continuum and Feiiopt modeling as well as subtraction. The 
specf it task was employed to make a reliable model of the Hfi 
spectral region. The simultaneous fitting of continuum and Feii 
should be less dependent on subjective evaluations by the ob- 
server We assume that the continuum underlying the Hy6 spectral 
region is a power law of variable slope. Continuum subtraction 
has severe limitations due to the small bandwidth of our spectra 
and to (unknown) internal reddening effects. However, specf it 
usually found a plausible continuum (see Fig.[T]i- Estimation of 
continuum subtraction uncertainty due to S/N involves choosing 
continuum fluxes at about -3cr (minimum) and +3cr (maximum) 
levels where <t is the standard deviation of the most likely con- 
tinuum choice. Uncertainties of continuum placement were em- 
pirically defined from the difference between extreme high/low 
continua and the most probable one derived by specf it. 

Feiiopt emission was modeled us ing a scaled and broad ened 
template as previously employed bv lMarziani et al.l ( 12003 ah . An 
important change was introduced to the Feii template emission 
underlying H/3 as described below ( 92.31 ). We stress again that 
a notable advantage of specfit is that the scaled and broad- 
ened template can be simultaneously fit over » 1000 A with 
a power-law continuum, the H/3 profile and the narrow lines 
(HjSnc, [Oiii]^^4959,5007, Fe vii15160, Feviil5177). The strong 
Fen blue blend (4450-4680 A: Feii/14570) was measured as an 
estimator of the Feiiopt strength. The Gaussian broadening fac- 
tor from the best fitting template yields an estimate of FWHM 
Fell. A careful estimate of minimum and maximum plausible 
broadening factors was made to derive +3cr uncertainties. The 
blue side of the spectrum including Feii/14570 is missing or only 
marginally covered in several sources. In these cases the best 
template fit was achieved for the red blend (Feii/i5130 in the 
range 5200-5600 A) and Feii/14570 was estimated assuming a 
fixed ratio between the red and blue blends. Fig. |2] shows the 
estimated Feii emission (green lines). 

The H/3 profile was modelled with: (la) a core Lorentzian 
component plus a weaker/broader Gaussian on the blue (usually) 
or on the red side of the core or, alternatively, (lb) the sum of two 
Gaussians one always very broad and showing non-negligible 
velocity shift relative to a narrower core component; (2) a narrow 
component (Hj6nc)- We do not attribute any physical meaning to 
the "blind" decompositions performed with specfit in the case 
of the double Gaussian fits for individual sources. The two types 
of fit provide the simplest empirical description of broad H/? pro- 
files in the majority of sources. The Hy6 profile was fit with a high 
order spline function (IRAF task sf it) in a few cases where the 



model fits showed significant deviations from the actual profile. 
Spline functions do not yield a model fit but only an empiri- 
cal fit that smooths noise and reproduces the main features and 
inflections in H/?. The spline fit approach was the one used in 
our previous work dMarziani et al l h003ah where sources with 
FWHM H/3 < and > 4000 km s ' favored Lorentzian and double 
Gaussian fits respectively. Systematic differences in flux and 
W measurements as well as in line profile parameters between 
the empirical and specfit approach occur since the extended 
wings of Gaussian and especially Lorentzian profiles lead to a 
lower continuum placement. The effect however is within the 
estimated measurement uncertainty. 

The [Oiii]/l/l4959,5007 lines were also simultaneously mod- 
elled with the following conditions strictly enforced during the 
fitting procedure: (1) the flux ratio between [Oiii]/I5007 and 
[Oiii]/l4959 ~3, (2) both lines should show identical profiles. It 
was often necessary to include a "semi-broad" and/or blueshifted 
component to model the [Oiii]/l/l4959,5007 profile since evi- 
dence for a strong blue asymmetry was observed in several 
cases. This approach yields a satisfactory reproduction of the 
observed [Oiii]/l/l4959,5007 profiles within our resolution and 
S/N limits. It is again an empirical method to properly estimate 
the total [Oiii]/l/l4959,5007 strength and describe the line profile 
shape. We do not make any assumptions about the actual nature 
of [Oiii]/l/l4959,5007 emission which will be the subject of a 
separate paper. 

2.3. The Feiiopt template 

Considerable attention has been given to theoretical and em- 
pirical estimation of Feiiopt emission in the spectral region 
of HyS, for a number of scientific (i.e., obtain a diagnostic 
of the BLR) and technical goals. Theoretical calculations of 
Fell emission, assuming the predominance of photoionization 
in the Broad Line Region have greatly improved in the last 
decade dSigut & PradhanL l2003t ISigut et al.L l2004t IVerner et al.L 
Il999h with a 37 1 -level model of the Fe^ ion pres ently included 
in th e photionization code CLOUDY v. 07.01 ( iFerland et al.L 
Il998h. A detailed analy sis of the Feii emission in I Zw 1 
( Veron-Cettv et al.l |2004|) revealed rich and complex emission 
from several line systems, each associated with a well-defined 
redshift and FWHM. In the case of sources with larger FWHM 
and data with intermediate resolution as wel l as average S/N an 
analysis as detailed as that made for I Zw 1 (IVeron-Cettv et all 
l2004l) cannot be done. The Feiiopt template approach works 
sat isfactorily for the he t erogeneous var i ety of type- 1 AGN 
(B oroson & Greenl 119921; iMarziani et all l2003at IZamfir et all 
2008^ with the probable exception of some outliers and ex- 
tremely strong Fell emitters. 

There is however a major problem associated with the 
empirical template: it is difficult to estimate the Feii emis- 
sion u nder the H/? line. Our em pirical template differs from 
that of iBoroson & Gree n (1992 ^ because o f lowe r Feii emis- 
sion under H/3 and from , Veron-Cettv et all (|2004|) for the op- 
posite reason. Th e new analysis of the I Zw I spectrum by 
IVeron-Cettv et al.l ( |2004|) . coupled with photoionization mod- 
els, suggests that emission under HyS should not be very strong. 
Therefore, the strength of the Feii lines underlying H/? has been 
computed with an appropriate photoionization model. In the 
context of photoionization the present understanding of the Feii 
emitting region suggests low ionization and high electron den- 
sity ( Verner et al., 1999; Bruhweiler & Verner, 2008)). We as- 
sumed an ionization parameter log F = -2.25, and electron den- 
sity [cm""*] log«e - 12.75. This model provides an intensity ra- 
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tio between Fen m29 and Hj8 consistent with that measured in I 
Zwl. These improvements in the template do not significantly 
affect our 7?Fen= F(Feii^4570)/F(H/?) ^ W(Feia4570)/W(Hy8) 
measurements. Even in the strongest Feii emitters (/?Feii~ 1) the 
change in H/J intensity is much lower than the associated un- 
certainty. The effect on H/? profile is to lower the Feiiopt sub- 
traction on the red side of Hy6 making it stronger. The change 
with respect to the previous template is negligible in the case of 
weak emission but may not be so if Feii is strong and relatively 
narrow (see 94.21). A blueshifted, broad component ascribed to 
a wind dSulentic et al.L 120021: iMarziani et al.L |2008|) is affected 
qualitatively by the new Feiiopt subtraction, although there are 
several examples of individual ISAAC sources (HE0248-3628, 
HE05 12-3329 in Paper II) that unambiguously confirm the ex- 
istence of such a component which is obvious in Lyg and es- 
pecia lly Civ/11549 for several low-z sources (IMarziani et al.L 
I2008h . 

3. Reduced Spectra and Line Measures 

3.1. Redshift Determination and Rest Frame Corrections 

Rest frame determination is important in quasar spectroscopy 
because both broad and narrow emission lines are known 
to show velocity displacements r elative to one another. We 
followed our standard procedure (IMarziani et al. I I2003ah and 
adopted the H/3 peak redshift [A = 4861.33 A) as rest frame 
for our new sources. The [Oiii]/l.i4959,5007 Unes [A 4958.9, 
5006.85 A) were only used when they gave the same value as 
HjS (see Col. (6) of Tab. [TJ. They sometimes show a blueshift 
relative to the H/ j peak and are designated as "blue o utliers" 
in extreme cases dZamanov et al .'. 2002; Marziani et al.. '2003bt 
lAoki et all 120051: iBorosonL [2005; Komossa et al.. .20081) . Blue 
outhers usually show FWHM Hj8< 2000 km s ' but caution 
should be used when adopting [Oin]/}5007 or Hy6 for rest frame 
determination (see also Hu et al., 2008a, who suggest a system- 
atic Feiiopt redshift with respect to the quasar rest frame). The 
adopted z estimates were used to de-redshift the spectra while 
the IRAF dopcor task applied a (1 + z)^ correction to convert 
observed specific fluxes into rest frame values. Fig.[T]shows the 
flux and wavelength calibrated spectra. 

3.2. Optical Feu Lines 

Fig.|2]shows the continuum subtracted spectra (left panels) with 
the adopted Feiiopt template model superimposed. Tab. |2] pro- 
vides flux and equivalent width of H/3, as well as flux and 
FWHM of Feii/14570. The data format is explained in the foot- 
notes of each table. Care should be applied in using the FWHM 
values for Feii/14570. These values were obtained through a 
minimization in the specfit procedure and therefore should at 
least be free of any subjective bias. Uncertainties were estimated 
by superimposing a sequence of templates with different FWHM 
values on the optimal Feii/14570 spectrum and estimating the 
±3cr confidence levels where FWHM Feii became unacceptably 
large or small. We find that FWHM(Feii/i4570) is systematically 
smaller than FWHM(HyS) and that the median value is a 4000 
km s"'. Uncertainties of individual Feii measurements are quite 
large allowing only statistical inferences to be made. 

3.3. Narrow Lines 

Table 13 reports the flux and equivalent widths values for HyGNC 
and [Oiii]/I5007. HjSnc is weak/marginal in most sources and 



undetected in 4 sources. The total flux values for [Oiii]/15007 
often include narrow and semi-broad components which are re- 
quired to adequately fit the [Oiii]/l/l4959,5007 line profiles using 
specfit. Differences between the UJSnc and [Om]/i/i4959,5007 
profiles indicate that crude intensity ratios [Oiii]/l5007/Hy6Nc 
may not be meaningful. The ratio [Oiii]/15007/Hj6nc has a clear 
meaning only if there is a redshift and FWHM consistency be- 
tween [Oiii]^^4959,5007 and H/3. 

3.4. n/3 Line Profiles 

Measurements of the broad HyS profile including FWHM and 
other important line parameters such as asymmetry index, kur- 
tosis and line centroid at various fractional intensities were ob- 
tained using a FORTRAN program developed for the purpose. 
These parameters are the same as defined in Paper I and are given 
in Tables|4|and|5] The right-hand panels of Fig. |2| emphasize the 
shape of H/3 for each source. Table[5]lists measurements of the 
Hy^BC centroid at various fractional intensities (in km s"'). Each 
line measurement in Tables |4| and |5] is followed in the next col- 
umn by its appropriate uncertainty at the 2cr confidence level. 
The uncertainties on line profile parameters were estimated by 
changing the fractional height by +0.05. All of these measures 
are affected by the compositeness of the broad H/3 profile which 
frequently shows unambiguous evidence for at least two distinct 
components: a broad component proper (Hp^c) and very broad 
component (HySvBC presumably from a very broad line region, 
VBLR). In the following, we will continue to refer to HyS as the 
total broad emission excluding H;Snc, and to HySBC as to the H/3 
emission after HySvBC (and of course also H/^nc) removal. 

4. Results 

4. 1 . Luminosity Effects 

H/i - Fig. [3] shows the distribution of H/3 profile measures as 
a function of bolometric luminosity where Lboi ~ 10 ■ AL,i 
{A = 5100 A) for the SDSS data. Bolometric luminosity was 
computed from z and for the ISAAC sample because of 
the spectra uncertain flux scale. Figure |3^ (upper left panel) 
shows the distribution of FWHM(Hy6) measures as a function 
of log Lboi where lo w z measures come from the SDSS sample 
(IZamfiretal.Ll2008h and ISAAC sources are indicated by larger 
filled circles. The only obvious, well defined FWHM trend found 
over ~5 decades of source luminosity involves an increase in the 
minimum FWHM H/3 with luminosity as previously reported 
in Paper I. Fig|3^ shows a systematic increase with log Lboi of 
mean/median FWHM H/3 from mean/median 3200/2550 km s"' 
to 5075/4900 km s"' which is likely driven by the minimum 
FWHM trend. No trend is observed, for example, among sources 
with FWHM > 4000 km s ' . The excess of sources with FWHM 
> 8000 km s"' in the range log Lboi = 45.5 - 46.5 simply re- 
flects Malmquist bias amplified by the increase of sources above 
log Lboi ~ 45.5 due to the onset of strong source evolution in the 
range z - 0.5 - 0.7. The correlation of FWHM H/3 with source 
luminosity is weak (c.f. Paper I) but statistically significant for 
the increased samples of this paper, with Pearsson's correlation 
coeflicient r * 0.26 (P £ IQ-^). 

Figure |3j' (upper right) shows the distribution of asymmetry 
index (A.I.) as a function of luminosity. An excess of red asym- 
metric profiles is seen at all luminosities but does not show up 
well in the mean or median values which lie in the range 0.03 - 
0.06 because of the strong SDSS source concentration between 
A.I. = + 0.05. The ISAAC sample shows the highest (reddest) 
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mean and median values with A.I. =0.11 and 0.09, respectively. 
If we focus on the sources with more extreme A.I. values then we 
find the largest difference between the "strips" -0.4 <A.I.< -0.2 
and 0.2 < A.I.< 0.4; red asymmetries outnumber blue ones by a 
factor of seven. Only 17 SDSS sources (zero ISAAC) show blue 
asymmetries greater than -0.2 while show red asymmetries 
greater than +0.2 (17 ISAAC). No luminosity correlation is de- 
tected with the red asymmetric excess visible at all luminosities. 

Figures [3}; (lower left) and[3}l (lower right) show the distri- 
butions of HySfic "peak" c(|) and "base" c(|) profile velocity dis- 
placements (line shifts), respectively, c(^) shows a strong source 
concentration around zero (unshifted profile base) but a signif- 
icant excess of redshifts at all luminosities. There is evidence 
for a trend with mean/median values increasing from 56/57 
km s"' in the lowest luminosity decade to 671/292 km s"' for the 
ISAAC sample. Thirteen sources (7 ISAAC) show base redshifts 
between 2000^000 km s"' while no sources show blueshifts in 
that range. Clearly the red excess in Figures [3j5 and [3]i are re- 
lated. If we distinguish between radio quiet (RQ) and radio loud 
(RL) sources (Zamfir et al., in preparation) we find the strik- 
ing absence of a RL source concentration near zero. Radio loud 
sources span the full range of observed red and blueshifts how- 
ever with a preference for redshifts. The concentration around 
zero is a pure RQ effect. 

The centroid c(|) shows only a small number of sources with 
shifts above +1000 km s ' (4 red and 7 blue and ISAAC). The 
approximately 2 for 1 preference for blueshifts is confirmed in 
the range c(|) - 600 -1000 km s ' with 8 red and 15 blueshifted 
sources. Mean values are slightly blue for lower luminosities (- 
20 to ^0 km s"' for log Lboi ~ 46) and become slightly red in 
the highest luminosity SDSS and ISAAC bins (+35 km s"'). The 
simplest accretion disk models predict redshifted profile bases 
and b lueshifted peaks ( Chen et a l.. 1989; Eracleous & Halpernl 
I2003h leaving open the possibility that an underlying disk signa- 
ture is present in all sources. Later discussion however does not 
support that conclusion for redshifted and red asymmetric pro- 
files. The absence of strong luminosity correlations in Figure [3] 
is confirmed by Pearson and Spearman correlation coefficients. 



[Oiii]/l/l4959,5007- Perhaps the most striking difference be- 
tween the SDSS and ISAAC samples involves the relative weak- 
ness of narrow line emission in ^ 80% of the latter. Table [3] 
shows only three ISA AC sources with W([Oi iiU5007) ^ 13 A. 
We previously found dSulentic et al.L l2000bh mean values for 
low z sources in the range W([Oiii]/l5007)= 23 - 27 A with stan- 
dard deviations for various subsamples in the range 10 - 30 A. 
Only one source in the ISAAC sample shows W([Oiii]/l5007) 
greater than the mean of the low z sample: HE 0109-3518 with 
W([Om]/15007) a! 38.5 A, an extremely large value especially 
considering the high luminosity of this z ~ 2.4 quasar This 
is the ISAAC source (Figure [1]) that most close ly resembles a 
low z quasar Twelve blue outliers identified by iMarziani et al.l 
(I2003bl) show W([Oiii]^5007) ^ 7.5A which is the same as the 
mean value for the sources listed in Table[3] There is a high frac- 
tion of sources in the ISAAC sample that may show a systematic 
blueshift with respect to Hj6nc, making several ISAAC sources 
similar to the blue outliers in terms of the [Oiii]/l/l4959,5007 
properties. The issue will be however discussed in a separate pa- 
per. 



4.2. 4DE1 Parameters 

A 4DE1 spectroscopic par ameter space previously defined 
dSulentic et al.Ll2000bll2007l) involves measures of FWHM(H/3), 
^F eii Civ/il549 shift, an d soft-X ray spectral slope (see also 
iBoroson & Greenlll992l) . We focus here on the optical plane of 
4DE1 where, as far as we can tell, none of the four principal pa- 
rameters show a correlation with optical or UV luminosity out to 
z ~ 0.8. Fig.|4]shows the SDSS and ISAAC source distributions 
in the 4DE1 optical plane (FWHM H/? vs. /?Feii)- Low-z quasars 
are distributed along a curved sequence with sources showing 
both large FWHM HyS and iJpeii apparently forbidden. The dis- 
tribution of SDSS sources matches closely the one found in the 
ATLAS sample. The zone of occupation extends from the upper 
left where FWHM(H^) is large, Feiiopt is least prominent and 
~50% of the sources are RL. At the opposite (lower right) end 
we find sources with the narrowest FWHM H/3 profiles, strong 
Feiiopt and ~1% radio loud sources. The centroid of the SDSS 
source distribution Ues near FWHM ^8= 3800 - 4000 km s ' 
and /?Feii= 0.3 - 0.4. 

We use the 4DE1 optical plane to see if the higher luminosity 
ISAAC sample shows any difference in source distribution that 
might hint at a luminosity trend. ISAAC sources are shown as 
large dots with error bars in Fig.|4] The largest distribution differ- 
ences involve: (a) the absence of ISAAC sources below FWHM 
~ 3000 km s ' and (b) an apparent displacement of almost half 
of the ISAAC sample towards larger values of FWHM H/3bc and 
/?Feii: on average by ^1 - 2000 km s"' in FWHM H/?bc and by 
0. 1-0.2 in /?Feii- Difference (a) involves the previously discussed 
increase in minimum FWHM H/?bc with luminosity and can be 
described as a zone of avoidance for high luminosity sources in 
the 4DE1 optical plane. Difference (b) requires more detailed 
study in the next section. 



4.3. 4DE1 Spectral Types and Composite Profiles 

Further analysis of 4DE1 and source luminosity trends are best 
accomplished using composite spectra. This avoids confusion 
introduced by details in individual spectra and optimizes S/N. 
Composite spectra in the context of the optical 4DE1 plane have 
special value because Feiiopt prominence is one of the key pa- 
rameters and FWHM H/3 is arguably the most reliable virial 
estimator of black hole mass. Fig. |4] is subdivided into bins 
of AFWHM(Hg) ^ 4000 km s ' and A Rp^n^ 0.5, following 
ISulentic et al](l2002l) . 

Spectral types (bins) Bl^, B4 and A3 each contain a sin- 
gle ISAAC source. A reexamination of the single B4 source 
HE1505+0212 suggests that Fen/14570 emission is unusually 
strong with no reasonable continuum adjustment capable of 
moving it further to the left than bin B3. This source is ei- 
ther a genuine outlier or an extreme Pop. A source (following 
the luminosity-dependent boundary between Pop. A and B dis- 
cussed in 34.5.1b . Many ISAAC sources are located in bin B2 
while SDSS sources are apparently rarer there and much more 
concentrated toward the lower left corner of the bin. R^eu is rel- 
atively insensitive to small changes in the adopted continuum 
leading us to conclude that the displacement of ISAAC sources 
towards higher /Jpeii cannot only be attributed to an effect of the 
data reduction. This supports the reality of distribution difference 
(b) discussed above. 



6 



Marziani et al.: in Intermediate z Quasars: III. Protile Analysis 



4.3.1. Decomposition and Fitting of ISAAC Spectra in tine 
Optical Plane of 4DE1 

Median composite ISAAC spectra were computed for sources 
in 4DE1 bins Al, A2, Bl, B2 (Fig. |5]l using calibrated and 
de-redshifted spectra normalized to specific flux of unity at 
5100 A. The first two composites show broader FWHM(Hy8) 
and the second two show stronger TJpeii compared to ATLAS 
jSulentic et al.Ll2002l) or SDSS composites for the same bins. 

Pop. A sources can usually be modelled with a single un- 
shifted Lorentzian component. This is also true for bin Al 
and A2 ISAAC sources. Pop. B sources usually require a dou- 
ble Gaussian (unshi fted broad Hj^Bc redshifted very broad 
HySvBC components: ISulentic et al.Ll2002h . Hy6 profiles in bin Bl 
and B2 sources are best fit with a double Gaussian model im- 
plying that both unshifted HySec and redshifted HySvBC compo- 
nents seen in low z quasars are also present in high luminosity 
sources (Fig.|6]l. The decomposition is not always unique and it 
can therefore be difficult to estimate the relative contributions of 
the two components. An important clue about the proper decom- 
posit ion comes from profi le variations in the source PG1416- 
129 (ISulentic et al.L |2000'3 see also 3C 206 in ICorbin & SmithI 
(I2OOOI) ) where the HySfiC almost disappeared following a con- 
tinuum decline while the H/3vbc showed less variation and 
was revealed much more clearly as very broad and redshifted. 
The highest S/N inflected spectra in the ISAAC sample (e.g. 
HEO 109-35 18) also clearly favor this decomposition which is 
adopted throughout our analysis. 



4.4. Luminosity Binned Spectra 

Comparing mean or median spectra for random samples of 
quasars observed in different luminosity ranges is unlikely to 
yield meaningful results. Any difference would be obscured by 
the intrinsic dispersion of spectral properties at any given ap- 
parent luminosity. As stresse d in much of our pre vious work 
(see especially the review by ISulentic et al.L l2000ah . and illus- 
trated again in Fig. |5] all quasar spectra are not self-similar 
Figure [3^ reenforces this point in one dimension by showing 
that FWHM H/3 spans a very large range from x 500 km s ' 
to R! 20000 km s"'. We suggest that a minimum requirement 
for physically meaningful comparisons among low-z samples 
involves distinguishing between 4DE1 Pop. A and B sources. 
More refined comparisons becomes possible by generating com- 
posite spectra for individual Pop. A and B bins: e.g. Al, A2, A3, 
. . . and Bl, B2, Bl^, etc. The population A - B concept was per- 
haps the most controversial part of our 4DE1 formalism where 
sources were arbitrarily divided into population A and B based 
on FWHM Hj6 profile measures above and below w 4000 km s ' 
respectively. However we find a multitude of differe nces when 
we distinguish sources this way (ISulentic et al.Ll2007l) . Most RL 
sources belong to Pop. B, and about 25-30% of RQ sources also 
belong to Pop. B, but the Pop. A/B distinction appears to be more 
effective than more traditional RQ vs. RL comparisons in terms 
of BLR structure. 

Fig. I2] shows median composite spectra in six luminosity 
intervals using 321 SDSS and 53 ISAAC sources. Composites 
were generated separately for Pop. A and B sources in decades 
of logLboi = 43 - 44 to 48 - 49. All ISAAC sources lie in 
the interval log Lboi - A1-A9 with a small number creating the 
highest Lboi bin (48-49). Pop A and B quasars show parallel 
trends with increasing luminosity: [Oiii]/l/l4959,5007 and H/^bc 
systematically decrease in equivalent width but HyS at the same 
time increasingly becomes redward asymmetric. The strong de- 



crease in narro w line strength is expected from st andard nebular 
theory (see e.g. lNetze3.ll990l with the caveats of ISulentic et al.l 
(I2OO7I) ). The Eigenvector 2 anti-correlation betwee n strength of 
Heii/14686 and luminosity (iBoroson & GreenifT992h is also seen. 
No obvious trend involving Feiiopt is seen which reflects the or- 
thogonality of luminosity with respect to Eigenvector 1 measures 
using principal component analysis (PCA). The ISAAC sample 
median spectra are most similar to the highest luminosity SDSS 
bin composite. The change in Hy6 towards high luminosity ap- 
pears to be concentrated in the strength of the red side of the line 
which was already anticipated as a red asymmetry and line base 
c{j) red shift in Figures |3]3 and [3};. 



4.4.1. Decomposition of Luminosity Binned Spectra 

The next logical step is to decompose and quantify individual 
lines and line components in different Pop. A and B luminos- 
ity bins. Earlier attempts revealed one of the most obvious dif- 
ferences among low redshift HjS profiles and p rovided one of 
the motivations for the population A-B concept ( Sulentic et al.L 



l2002l:lMarziani et al.Ll2003bl:ISulentic et al.Ll2007l) . Fig.l8lshows 
Pop A and B median spectra binned in luminosity after contin- 
uum subtraction. Tab. |6]provides line equivalent width measures 
and FWHM Feii values for various spectral bins discussed in the 
text. Tab. [T] tabulates line profile measures and Tab. |8] provides 
H/3 centroid measures and estimated uncertainties for the same 
bins. Formats appear as footnotes below each table. The first 
seven lines present various 4DE1 bin values for ISAAC sources 
while the remainder provide luminosity bin values following the 
designations in Fig.|4| Measures in these tables reflect the results 
of simultaneous fits to all lines and components in each median 
composite spectrum. Attempts at a sequential line-by-line de- 
composition yielded very similar results. Numbers quoted in the 
discussion of Figure |3] are averages or medians based on indi- 
vidual source measures while Tables |6]|7]|8] give values derived 
from composite spectra. There are no significant discrepancies 
between measures estimated in these two ways. Source codes 
for Lboi bins use the first two digits of the Lboi bin followed by 
A or B for the population (e.g. values for line 43A correspond to 
measures of the Pop. A luminosity composite for log Lboi = 43 
- 44). Narrow line changes with luminosity will be discussed in 
a separate paper so that we can concentrate on the properties of 
Hj6bc and HySvec- FWHM Hj6 in Pop. A sources shows a system- 
atic increase with Lboi from 2000 km s"' to 4100 km s"' across 
the six luminosity bins likely driven by the increase in minimum 
FWHM Hj6 with log Lboi shown in Figure [3^. 

Pop. B shows no clear FWHM trend with log Lboi but we find 
that W(Hj6) decreases systematically with increasing luminos- 
ity. Pop. A does not show an obvious trend although W(Hy6) is 
smallest in the most luminous Pop. A bins. The Pop. B "Baldwin 
Effect" shows Wilifi) decreasing from 143 A to 82 A with in- 
creasing log Lboi. Pop. B sources show 6-7 times larger asym- 
metry indices at the highest source luminosities as well as evi- 
dence for a redshift near the base (c(|)) of the line. The mean 
velocity redshift increases by 9 times from the lowest to highest 
luminosity bin. Fig.|9|shows the increase with source luminosity 
of the fraction of c(i) redshifts for Pop. B sources. The excess 
flux in the red wing of Hj8 is revealed in Fig.|8]decompositions to 
be a real H;3 flux excess and not du e to Feiiopt contamination (see 
also lNetzer & Trakhtenbroll l2007l who found a similar result). 

Conversely Feiiopt subtraction reveals that shows little 
evidence for profile shifts and asymmetries in Pop. A sources. 
The apparent Pop. A red asymmetry appears to increase with 



Marziani et al.: Hj8 in Intermediate z Quasars: III. Protile Analysis 



7 



luminosity in Figure |2l however Figure [8] suggests that it is 
mainly due to Fenopt and [Om]/l/l4959,5007 contamination on 
the red side of H/?. Table shows that W(Feii/i4570) does not 
increase with luminosity but H^Sbc decreases somewhat mak- 
ing red Feiiopt contamination appear to grow in strength, even 
if some redward, very broad emission cannot be excluded at a 
low flux level in the highest luminosity bins. 

The quantification in Tables |6] |2l [H and the two compo- 
nents required to model Pop. B sources, ree nforce past evidence 
(ICorbinl[T99l ISulentic et all l2000d l2002h for the existence of 
a redshifted very broad line component in HjS. The HySvBC is 
the most obvious profile difference between Pop. A and B, as 
well as between high and low luminosity sources in Pop. B. As 
noted earlier both RL and RQ sources are found in Pop. B and 
both Pop. B RL and Pop. B RQ sources show the HjSvbc- This 
provides a strong spectroscopic argument in favor of the impor- 
tance of the Pop. A-B concept compared to comparisons on the 
basis of radio loudness. The effect of the VBLR component is to 
produce a broader, more red asymmetric and redshifted H/J pro- 
file with the amplitude of these effects dependent on the relative 
strengths of the BLR and VBLR components. HE 0109-3518 in 
Fig. [U is an example of a source with a weak/moderate VBLR 
component while HE 1039-0724 shows one of the strongest 
VBLR components in the ISAAC sample. Tab.|7][8l Fig.|9l and 
the profile decompositions in Figure [8] show evidence for a uni- 
form increase of H/?vbc strength with luminosity. 

4.5. Luminosity Dependence of Spectral Types 

The left-hand panel of Fig. [TO] shows the bolometric luminos- 
ity distributions for ISAAC sources showing H/3 profile best fit 
with single Lorentzian and double Gaussian models. Ambiguous 
sources are excluded. A K-S test confirms that the two source 
luminosity distributions are not significantly different. The right 
panel shows the distribution of profile fits in an FWHM li/3 vs. 
bolometric luminosity plane with Lorentzian, double-Gaussian 
and ambiguous sources denoted L, G and U respectively. No 
strong luminosity preference is seen for L or G sources. The 
difference between FWHM distributions for L and G profiles is 
highly significant with a K-S test suggesting that P ^ 10"** for 
sources with different profile shapes to be are drawn from the 
same parent population. This can be taken as an independent 
confirmation of the Pop. A-B dichotomy. 

4.5. L Spectral Types: A Refinement 

The Pop. A-B FWHM(Hj6bc) = 4000 km s ' boundary was 
adopted using a low z and low luminosity sample. There is no 
strong evidence that it has a direct physical meaning even if 
it is better motivated empirically than the limit at FWHM(Hj6) 
= 2000 km s"' used to separate "Narrow Line" Seyfert Is 
(NLSyls) from the rest of broad fine AGN. 

Fig[TT]presents a reprise of Figure[3h with a much expanded 
low redshift sample including: (a) SDSS DR5 quasars brighter 
than g=17.0 within z=0.7, (b) FIRST radio dete cted SPSS 
DR5 quasars between g=17.0-17.5 within z=0.7 (Za mfir et al.L 
l2008h and (c) the AT LAS sample of low redshift quasars 
dMarziani et al.Ll2"003ah . This provides a somewhat better map- 
ping of the low z source distribution in the FWHM - log Lboi 
plane. The SDSS quasar sample does not include sources with 
FWHM HyS<1000 km s ' which were classified in SDSS as 
galaxies. The narrowest of these "broad line" sources can be con- 
firmed as Type 1 AGN (NLSyl) by the presence of Feii/14570 



emission in their spectra. Among the 4 00-t- highest S/N SD SS 
quasars that we have studied in detail (IZ amfir et all l2008h all 
but five show detectable Feiiopt emission suggesting that the pres- 
ence of optical Feiiopt emis sion is an ubiquitou s property of Type 
1 AGN. A recent survey dZhou et al.L |2006|) of NLSyls using 
SDSS DR3 provides a representative census of these sources. 
They found 41 sources (z $ 0.7) brighter than g = 17.5 be- 
tween FWHM H/3= 450 - 2000 km s ' and these are indicated in 
Figure[TT]with filled triangles. The old Pop. A-B and NLSyl -rest 
of Type-1 AGN "boundaries" at 4000 and 2000 km s"', respec- 
tively are indicated. This figure provides us with a better defined 
lower edge to the FWHM Hy6 distri bution. We superpo se curves 
corresponding to exponents a=0.52 dBentz et"alll2006l ) and 0.67 
(Paper II) in the Kaspi relation which will be considered in the 
discussion. Fig[TT]also shows a revised Pop. A-B boundary fol- 
lowing the form of an a = 0.67 Kaspi relation as done in Paper 
I. 

Figures[TO]and[TT]suggest a modified 4DE1 Pop. A/B bound- 
ary is required for high luminosity quasars. We define modi- 
fied 4DE1 spectral types MAI, MA2, MBl, MB2, etc. where 
the boundary between the MA and MB Populations is now 
luminosity-dependent and given by FWHMab ~ 3500 + 625 ■ 
^boi^^ (a - 0.67; analogous to Paper I). Composite line li/3 pro- 
files for the revised 4DE1 bins are similar to those obtained with 
the fixed FWHM definition (Fig. |5] and Fig.|6l) and are therefore 
not shown. 

In order to decide whether the fixed (4000 km s"') or lumi- 
nosity dependent boundary is more significant we defined three 
samples: (1) the old Pop. A, (2) an intermediate population be- 
tween 4000 km s ' and new luminosity-dependent hmit and (3) 
a modified Pop. B MB = B \ M. We can also consider a mod- 
ified Pop. A, MA s A U M. The H/3 spectral region before and 
after continuum subtraction for the median spectra of A, M and 
MB are shown in Fig. [12] and Fig. [13] respectively. The M me- 
dian most closely resembles the Pop. A median rather than MB 
(or B). Fig[T2] and [13] show that the M HyS spectral region re- 
sembles the one of Pop. A, and that the H/3 fit is well fit by a 
Lorentzian function, as it is the case for Pop. A sources. This 
suggests that the luminosity dependent limit FWHMab (i') might 
be more appropriate. The distribution of G and L sources in the 
FWHM(HySBc) vs. luminosity plane (Fig.[TQ]i also supports this 
result. Although individual fits should be treated with caution we 
see that the shapes are predominantly Lorentzian in the Pop. A 
and M regions while they are Gaussian in the MB area of the 
plane. 

4.6. Black Hole Mass and Eddington Ratio 

FWHM(HjS) has been used as a virial estimator of Mbh 
dVestergaard & Petersonl |2006|) at low redshift. The problem is 
what to use above redshift z ~ 0.9. The use of CiV/11549 as an Hp 
surrogate involves problems that seriously reduce the reliabil- 
ity of high z Mbh estimates dSulentic et al.L l2007t iNetzer et all 
|2007|) . We have chosen instead to follow H/3 out to the highest 
possible redshifts via infrared spectroscopy. We use the updated 
relationships linking FWHM(Hy6Bc), specific luminosity at 5100 
A and BLR radial distance Rblr from IVestergaard & Petersonl 
d2006l their Eq. (5)). We do not enter into the caveats of 
black hole mass and Eddington ratio determinations (see e.g . 
Collin et al. ', '2006!; ISulentic et all l2006t iMarconi et all 120081: 
Shen et al., 2008) remarking only that Mbh values are uncer- 
tain by a factor 2-3 at Icr confidence level if derived from 
single-epoch observations. It is also possible to infer bolomet- 
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ric luminosity from the specific continuum flux assuming a con- 
stant bolometric correction. One can therefore also estimate 
the luminosity-t o-mass ratio (Lboi/MBH)- This approach is very 
crude (see e.g. iHopkins et atl 120071; iKellv et all |2008) but a 
bolometric correction is still relatively stable across RQ and 
most RL sources with the obvious exception of core-dominated 
(beamed) sources. 

The distribution of ISAAC quasars with L- and G-type fits as 
a function of Mbh and Lboi / ^^Edd are shown in the left- and right- 
hand panels of Fig. [T4l respectively. Fig. [TS] shows all ISAAC 
sources (including U HySec profiles) in the Lboi/i'Edd vs. Mbh 
plane. We derive quite large masses in the range log Mbh=8.5- 
10.0 with a range of log Lboi/i'Edd= ^1-0 to 0.0. L (assumed 
Pop. A) and G (assumed Pop. B) profiles favor lower/higher log 
Mbh and higher/lower log Lboi/^^Edd val ues, respectively, as pre - 
viously found for our ATLAS sample jMarziani et al. I l2003bl) . 
The difference in the distribution is significant for both Mbh and 
^boi/i'Edd, with a low probability that sources with L- and G-type 
profiles come from the same parent distribution of Lboi/i'Edd and 
Mbh: K-S tests indicate a probability % However, there is 
a range of Lboi/i'Edd(-0.5-f -0.3) where L- and G- type sources 
are almost equally represented. The absence of a sharp Lboi/^Edd 
boundary can be understood in terms of error effects as well as of 
systemic biases. Typical uncertainty in logLboi/^^Edd is ~ 0.1, if 
one assumes that the bolometric luminosity is uncertain by 50%, 
and that the FWHM(H;6) has typical eiTor of 10%. However, this 
estimate neglects the scatter in the size-luminosity relationship 
for the BLR. In addition, orientation effec ts - which are estab- 
lished among radio loud quasars (see e.g.. lSulentic et all 120031: 
iRokakiet all 120031) - may well raise the Lboi/i'Edd of several 
sources of Pop. B. They are not considered in the Lboi/i^Edd and 
Mbh estimates presented here. Moreover, the magnitude limit of 
our ISAAC sample makes possible that we can detect only Pop. 
B sources with the largest Lboi/^Edd values (cf. Fig. 7 of Paper 
II). 

Fig- dH shows the dependence of c(|) and A.I. on Mbh and 
^boi/i'Edd- Significant correlations between Mbh and both A.I. 
and c(|) are found while no significant correlations are found 
involving Lboi/i'Edd- The A.I. vs Mbh correlation is marginally 
significant (P a; 4 ■ 10"^) if Pop. B sources are considered alone; 
in this case c(|) and Mbh are not significantly correlated. This is 
perhaps not surprising considering the small number of sources, 
the typical errors, the Mbh "segregation" between Pop. A and 
B. Median spectra in bins of Mbh and Lboi/i'Edd point toward 
the origin of the correlation with Mbh (Fig. [TtI i. The bin la- 
beled as log Mbh< 9.5 actually covers masses in the range 
8.5 ^ log Mbh< 9.5, and for logLboi/i'Edd^ - 0.5 corresponds 
to the b in 3.9 ^ log Lboi/ Mbh^ 4.4 in Figure 9 of lMarziani et al.l 
(I2003bh . The remaining three bins represent the extension to 
higher masses and Lboi/i'Edd that were not sampled at z ^ 0.7. 
We see a trend, as observed at low-z, that for fixed Lboi /^Edd the 
redward component is stronger for larger Mbh- This result is dis- 
cussed in 35.41 

5. Discussion 

5. 1 . Luminosity Trends 

One of the motivations for the ISAAC survey was to search 
for trends, or even correlations, between emission line proper- 
ties and source luminosity. Another motivation was to compare 
high and low luminosity source properties in the 4DE1 c ontext. 
PCA of a bright quasar sample dBoroson &GreenL[T992h found 
source luminosity to be orthogonal to all line measures except 



W(Hen/i4686). An [Oiii]/I5007 meas ure was included among 
the original Eigenvector 1 parameters ( iBoroson & Greenl Il992h 
showing a large scatter likely due to the strong Lboi dependence 
that we find with ISAAC. We did not adopt an [OIII] measure 
in our 4DE1 formalism because we pre ferred to avoid narrow 
line considerations as much as possible dSulentic et al.Ll2000bl) . 
Subsequent 4DE1 studies at low z (e.g.. Paper I and references 
found no luminosity dependence on broad line LIL measures. 
What was thought to be the most significant (inverse) corre- 
lation between lumi nosity and broad lin e HIL measures (es- 
peciafly Civ^l549) ([Baldwin et al.Lll978l) is now known to be 
present in quasar samples with small redshift/luminosity dis- 
persion dXu et all [20()8J and is lik ely related to intrins ic source 
properties such as Eddington ratio dBachev et alll2004h . 

Only one luminosity trend was found over the 5-6 luminos- 
ity decades covered in this study (without Pop. A/B discrimina- 
tion). Profile parameter distributions (Figure [3) are essentially 
the same at all luminosities. The one exception involves an in- 
crease in minimum FWHM HJB with log Lboi visible in Figure 
[3^ (see also Paper I). Low luminosity (log Lboi - 43.0) sources 
show FWHM HjSbc values as low as 500-750 km s"' while the 
minimum value at log Lboi = 48 appears to be FWHM Hfi * 3000 
km s ' . The lower FWHM H/3 limit is expected if broad Balmer 
line emission arises in a virialized clo ud distribution (or accre- 
tion disk) that obeys a Kaspi (see e.g. iKaspi et all l2005h rela- 
tion. We can think of the minimum FWHM-log Lboi trend as the 
boundary for sources radiating sub-Eddington. NLSyls defined 
(at low z) as sources with FWHM H/3 <2000 km s"' are absent 
from the ISAAC sample. 

It is unclear if a Kaspi relation {Rblr L^^) is valid above 
redshifts of a few tenths. There has been much discussion about 
the correct exponen t in the Kaspi relation with values between a 
= 0.5 - 0.7 favored (iKaspi et al.ll2000ll2005LlBentz et alll2006l) . 
Figure nn provides an empirical representation of the minimum 
FWHM Hy6 trend. Note that the expected trend has been quadrat- 
ically combined with the instrumental resolution (assumed to be 
200 km s ' at FWHM, appropriat e for the SPSS). S uperposed 
curves corresponding to a = 0.52 ( iBentz et all l2006l) . and 0.67 
(Kaspi et al., 2005, Paper II) show that discrimination between 
these values becomes possible above log Lboi ~ 47. A first im- 
pression is that a - 0.52 is disfavored by the ISAAC data unless a 
considerable number of these sources are super-Eddington radi- 
ators. Thus values in the range a = 0.6-0.7 are favored. However 
a complicating factor involves the dependence of FWHM Hfi on 
source orientation which is especially important if one accepts 
the paradigm of a line emitting accretion disk. Some low redshift 
sources with small FWHM H/3 (extreme NLS yl sources) show 
rapid high amplitude X-ray variability (e. g. iBrinkmann et all 
I2003h or blue shifted narrow lines (blue outliers) which have 
been a rgued to signa l near f ace-on orientation. In a previous 
paper dMarziani et all l2003bl) we attempted inclination correc- 
tions for blue outliers assuming a sin / dependence of FWHM. 
Using a tentative correction that increases black hole mass by 
AMbh~ 0.4 we were able to move many apparently superEd- 
dington NLSyl radiators (sources with log Lboi ~ 46 that lie be- 
low the a = 0.5 curve in FigurefTTT) below Lboi/LEdd= 1- It seems 
unlikely that more than a very few ISAAC sources involve pre- 
ferred orientation (5 are nominally super-Eddington) given the 
distribution of observed FWHM H/3 values even if optical lumi- 
nosity is orientation dependent. It is also possible that orientation 
becomes increasingly less important in higher Lboi sources be- 
cause of radiative instabiliy in the accretion disk (jBlaes, 2007]). 
A value as large as a = 0.67 would imply that all quasars radi- 
ate well below the Eddington limit. We conclude that there is no 
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strong evidence for a population of super-Eddington radiators in 
our low z and ISAAC samples. 



5.2. 4DE1 Populations A and B 

This study suggests that, in order to isolate luminosity effects 
and other fundamental trends in quasar samples, one should 
work in the context of the 4DE1 formalism and take advantage 
of the Population A-B concept. The empirical evidence contin- 
ues to grow in favor of two distinct populations of quasars. We 
recently summarized the m ultitude of ernpirical differ ences be- 
tween the two populations jSulentic et al.L l2007ll2008l) . If lumi- 
nosity trends exist for only one of the populations then the lack 
of a trend in the other might obscure the effect in a mixed sample 
such as the one displayed in Figure [TT] The Pop A-B bound- 
ary was empirically defined, after noting a discontinuity between 
the sources whose FWHM(Hy8) lies in the range 2000 - 4000 
km s 'and those of Pop. B (also because almost all RL sources 
are pop. B). The similarity in emission line parameters between 
those sources and the NLSy Is then suggested a single population 
for the entire 0-4000 km s"' range. In analogy with the minimum 
FWHM Hfi trend which corresponds to a fixed Eddington ratio 
of 1 we can expect that the limit between Pop. A and B will also 
be luminosity dependent (Paper I). 

We have previously suggested that the Pop. A-B FWHM 
HjSbc boundary might reflect a critical accretion rate/Eddington 
ratio {L\,oi/LEdd~ 0.2 + 0.1). We note that a number of recent 
studies have indepe ndently arrived at a similar "critical" FWHM 



or Lboi/LEdd value ([Marziani et all l20q3bl: ICollin et al.L 120061: 



iBonning et all 120071: iKeflv et al.L 120081: iHu et all l2o68bl) . Such 



a critical value might signal a significant change in disk struc- 
ture/kinematics that is reflected in systematically larger FWHM 
Hj6, the demise of a CiV/11549 wind as well as weaker soft 
X-ray and Fen emission. It is interesting to point out that the 
thin accretion disk approximation is expected to break down at 
Lboi/LEdd~ 0-3 with t he development of an inner thick structure 
dHubenv et aUl2000b . This kind of change might well affect the 
spectroscopic signature of the BLR. 

It is probably naive to expect a luminosity dependence as 
strong as the one predicted by a fixed Lboi/i'Edd (Fig. fTTTi. The 
Pop. A/B boundary is found to be close to 4000 km s"', fol- 
lowing four lines of evidence: (1) the median spectrum for Pop. 
A in the luminosity bin 43$ logLboi $44 shows that the Hy6 
is Lorentzian; (2) between constant FWHM(Hy6)=4000 km s"' 
and the Lboi/i'Edd= 0.15 line, and for 45 $ logLboi $ 47, 
most ATLAS sources show typical Pop. A properties. A minor- 
ity of core-dominated radio-loud sources are located there since 
they are expected to be observed almost pole-on, at minimum 
FWHM. If 44 $ logLboi $ 45 ATLAS sources are mainly lo- 
cated within -500 km s"' of FWHM = 4000 km s"', their H/3 
profiles are similar to but not exactly the same as the ones of 
Pop. A. In the present paper we consider the luminosity depen- 
dence of the boundary to follow Paper I assuming a x 0.67. It is 
assumed constant corresponding to FWHM~4000 km s"' below 
logL =s 47. 



5.3. An HjSbc "Baldwin Effect" is Present only in Population B 

The Pop. A/B concept provides a useful alternate way to evalu- 
ate luminosity properties. When we discriminate by source pop- 
ulation we immediately find possible luminosity trends between 
FWHM and W(Hfi) for Pop. A and Pop. B respectively. The for- 
mer trend sees a mean FWHM H/? increase with luminosity from 



a; 2000 km s ' up to » 4000 km s This is likely driven by 
the minimum FWHM H/3 trend already discussed as well as by 
Malmquist bias since all of the highest luminosity sources lie 
at higher redshifts and are preferentially selected in a flux lim- 
ited sample. Assuming an underlying correlation between black 
hole mass and source luminosity we expect systematically larger 
black hole masses for the ISAAC sources as is observed. 

Pop. B sources show a decrease in W(Hy6) by almost a factor 
of two over our source luminosity range. This Hfj Baldwin eff'ect 
is difficult to interpret because HyS is a composite feature (H/Sbc 
+ HjSvBc)- Tab|6] presents an attempt at quantitatively separat- 
ing BLR and VBLR components for the six Pop. B luminosity 
binned composite spectra. The same can be done for Pop. A but 
the BLR component always dominates H/? in all luminosity bins. 
The BLR component in Pop. B ranges from 60% of H/S at low lu- 
minosity to 25% at log Lboi = 47-48. Correction for the VBLR 
component yields a "super" Baldwin effect where W(Hy6Bc) de- 
creases by a factor of four over our luminosity range. At the 
same time we find considerable scatter but no obvious trend in 
the strength of Fen emission as measured in the Feii/14570 blend. 
The highest luminosity Pop. B sources with smallest W(H/6bc) 
show an apparently more prominent HjSvbc while its contri- 
bution in low luminosity sources seems modest because HJSbc 
dominates. 

The simplest interpretation of these results involves changes 
in physical conditions of the BLR emitting region due to: (1) 
source luminosity and (2) the Pop. A-B dichotomy. As alluded 
to above the first may involve changes related to radiation pres- 
sure dominance in the accretion disk at higher source luminos- 
ity and the second changes related to a critical Eddington ra- 
tio. The tendency for FWHM(Feii) to follow FWHM HjSbc in 
BLR-dominated Pop. A sources can be taken as evidence for 
an Fell - HJB B LR commonality. The inferred lower density of 
the VBLR (e.g., ISulentic et an.f2000c[) in (less BLR dominated) 
Pop. B so urces disfavors an Fe ii/VBLR emission region com- 
monality dMarziani et al.L l2008h . W(Hj6vbc) is found to be ap- 
proximately constant in the six luminosity bins suggesting that 
it arises largely in optically thick gas photoionized by the cen- 
tral continuum. The largest challenge is to explain the four- 
fold decrease in HjSbc strength observed only Pop. B sources. 
The simplest interpretation involves a change in Pop. B phys- 
ical conditions that suppresses HjSbc- Since HySBC and Feii are 
thought to arise in the same region this might point toward colli- 
sional quenchi ng of H/3 BLR emission as originally suggested by 
lGaskel]| ( IT988l) . Another, even simpler possibility involves differ- 
ent sized annuli for Feii and H/Jbc emitting regions. If we correct 
W(H/3) measures for HjSvbc and recalculate /Jpeii parameters we 
find that many/most Pop. B sources will move towards higher 
Rpsii values a:l in the 4DE1 optical plane. The similar /?Feii ra- 
tios observed in Pop. A luminosity bins (using H/Sbc only) sug- 
gest also a similarity in physical conditions. In other words, we 
would have a similar low-ionization BLR in both Pop. A and B 
sources plus a second VBLR only in Pop. B sources. 



5.4. Black Hole Mass and Eddington Ratio 

Tab. |6] presents FWHM H/^bc estimates for Pop. B luminosity 
bins. The H/3bc of Tab.|6]is related to the "reverberating" compo- 
nent almost equivalent to the entire H/3 profile in Pop. A sources. 
We note that the decomposition H^Sbc- HySvBC of Pop. B sources 
assumes two Gaussians while the intrinsic shapes of H^bc and 
H;6vBC, especially in the radial velocity range where they overlap 
is, unknown. The origin of the redward asymmetry often seen in 
the Hy6 profile is also unclear. It seems too large to be due to grav- 
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itational redshift associated with line emission from increasingly 
smaller distances from the central black hole (under the assump- 
tion that the motion is predominantly virial). A likely possibility 
is non-virial motion as sociated with gas infall toward the cen- 
tral b lack hole (see e.g. iHu et al.L l2008at lOaskell & Goosmamil 
I2OO8I who provide renewed support for the hypothesis of in- 
falling emitting regions), but alternati ves like Compton scatter- 
ing have also been considered (e.g. ,Kallman & Krolik, J986i). 
However, if the VBLR is mainly made of optically thick gas the 
radial velocity range underlying the H/3bc is especially ambigu- 
ous. Note also that the width reestimatio n is not the same as 
using the 2"* ^ moment of the li/3 profile: cr dPeterson et al.Ll2004t 
ICollinet all 12006) may have little meaning for a composite line 
since it is measured on the whole H/? profile. The reestimation 
leads to a more extreme reduction than suggested in Paper II; 
actually, all VBLR corrected FWHM values decrease to within 
1000 km S-' of the Pop. A-B boundary (4-5000 km s"'). In fact 
the highest luminosity bins show FWHM(Hy6Bc) within the re- 
vised Pop. A range as defined in Figure [TOl 

An important implication of this correction involves its effect 
on Mbh and Ledd estimation for Pop. B sources. Pop. B sources 
would essentially cease to exist as a separate population in the 
context of Mbh estimation. All Mbh values for Pop B would de- 
crease from the previous range (Paper 11) log Mbh= 7.5 - 9.5 to 
7.0 - 9.0 which is mo re similar to the Pop. A range but still con- 
sistent with the view dSulentic et a l.. 2000a) that Pop. B quasars 
are radiating at lower Lboi/^^Edd, as further discussed below. 

Considering the uncertain and somewhat speculative nature 
of the corrections applied to Pop. B sources, FWHM H/? in Pop. 
A sources appears to provide more direct and reliable black 
hole mass estimates yielding a range in log Mbh = 6.0 - 8.5 
at low z and 9.0 - 10.0 at z ^1.0 (see figures in Paper II). 
The major source of uncertainty involves the inclination cor- 
rection which depends on the BLR geometry. A weak VBLR 
component and/or blue asymmetry are present in an uncertain 
number of sources which might lead to slight overestimates of 
Mbh when assuming that the virial line component i s a symmet- 
ric Lo rentzian. Gaussian fits to Pop. A profiles (e.g. lShen et al.L 
l2008h lead to significant overestimates of Mbh- Median values 
from our luminosity composite profiles (assuming mean FWHM 
Hj6 and Lboi for sources in that bin) show a clear luminosity 
correlation ranging from log Mbh~ 6-10 for the most lumi- 
nous ISAAC sources. Estimates corrected for Hy6 asymmetries 
are slightly lower but only for the SDSS bins. 

The Lboi/i'Edd values in Table |6] increase from 0.02 to ^ 
0.6 in the highest luminosity bins of Pop. B. We remark that 
the Lboi/i'Edd value for the most luminous bin (^ 1) is espe- 
cially uncertain because of the limited number of sources and 
of the large correction ampHtude (FWHM(Hj8)a!9300 km s 
FWM(Hj6bc)~4300 km s"'). Lboi/^Edd values obtained using 
FWHM (Hj6) without removing HjSvbc (Table |7]) range from x 
0.02 to 0.27. Similar values are obtained if FWHM(Feii/l4570) 
is used as a virial estimator to compute Mbh- The corresponding 
Lboi/^Edd estimates for Pop. A luminosity composites also show 
an increase ranging from Lboi/i'Edd ~ 0.2 (lowest luminosity) to 
1.3 in the highest luminosity bin. Ranges for Pop A from Paper 
II were Lboi/^Edd = 0.05 - 1.0 at low z and 0.5 - 1.0 at z > 1.0 
assuming that sources above Lboi/i'Edd>l-0 are the result of face- 
on orientation. The luminosity trend is expected for both Pop. A 
and B since we are sampling preferentially higher Lboi/i-Edd radi- 
ators with increasing Lboi (as shown in Paper II). Without HjSvbc 
correction the difference between the Lboi/i-Edd of Pop. A and B 
is A log Lboi/i'Edd* 0.6 for the SDSS luminosity bin, decreasing 
to «; 0.3 for the ISAAC luminosity bins (47 and 48). If we con- 



sider the width of H^Sbc only, A log Lboi/i'Edd^ 0.5, even if this 
difference becomes rather small in the ISAAC bins, probably 
because of the extreme amplitude of the correction, and of the 
caveats listed at the beginning of this section. At any rate, from 
the SDSS sample we can conclude that Pop. B appears to be a 
population of lower Eddington ratio than Pop. A. Pop. B is likely 
to be an older and more evolved quasar population than Pop A, 
since black holes seem to reside in h osts with large sp heroidal 
components, often elliptical galaxies ( iWoo et al.Ll2005l) . On the 
contrary the morphology of NLSy 1 hosts (i.e., in many ways the 
"extreme Pop. A" sources) in the local Univers e often involves 
high surface brig htness star-forming galaxies dKrongold et al.L 
I2OOII: lOhta et"an. l2007i) . They are often barred and/or remark- 
ably perturbed suggesting that they may be young systems sus- 
tained by a large flow of matter toward thei r central black holes 
(see e.g. the review bv lSulentic etaill2008h . 

5.5. Occupation of Spectral Type B2 

ISAAC sources are found in the B2 bin of the 4DE1 optical plane 
(FWHM= 4-8000 km s-' and /?Feii= 0.5 - 1.0. The low z ATLAS 
sample included only 1 1 sources (~5%) in that bin but the more 
complete SDSS sample involves 40 B2 sources (~10%). Bin B2 
may involve sources with the largest black hole masses (several 
times 10^ M©). The largest masses are expe cted to fall in the 
region of bin B2 jZamanov & Marzianii |2002|) but such massive 
black holes are not observed at low z (see, for example the Mbh 
distribution as a function of z in Paper II). In the ISAAC sample 
we derive Mbh values as high as log Mrh ~ 9.7 tSulentic et afl 
2006). According to the grid in Fig. 2 of IZamanov & Marzianii 
(120021) the B2 bin is expected to be occupied by sources with 
Mbh ~ 5 X 10'' Mq. The expected Lboi /Mbh ratio is modest with 
Lboi/LEdd~ 0.16. B2 sources therefore expected to share physical 
properties that are thought to be typical of Pop. B. The HyS line 
profile is well described by a double Gaussian decomposition 
also if we restrict our attention to revised bin MB2. This result 
is consistent with the conclusion that most sources in bin B2 are 
similar to those in Bl although with stronger /?Feii- The H/3 line 
profiles imply that both types may fall on the Pop. B side of a 
critical Lboi/^Edd that separates Pop. A and Pop. -B boundary. 



6. Conclusions 

We presented VLT-ISAAC spectroscopic observations for 30 
intermediate z and high-luminosity HE quasars. Combined 
with previous data we have a sample of 53 objects (z = 
0.9 - 3.0) whi c h we compare with two l arge low-z samples 
(iMarziani et al.L l2003ai: IZamfir et all |2008|) . We find few cor- 
relations/trends between broad emission Une properties and 
source luminosity. We previously proposed an empirical limit 
at FWHM(Hy6) x 4000 km s"' to distinguish between two popu- 
lations A-B of quasars with very different spectroscopic proper- 
ties. Sources with FWHM(Hy6) < 4000 km s ' show Lorentzian 
Hfj profiles while those above this limit are best fit with dou- 
ble Gaussian models. This low z/ low luminosity dichotomy is 
also found for the high luminosity ISAAC sources. Elementary 
Lboi/^Edd computations suggest that the dichotomy may be ex- 
plained if a critical Lboi/Lfidd is associated with a BLR structure 
change, but the issue of the Pop.A/B boundary deserves further 
scrutiny. Once the population A/B dichotomy is taken into ac- 
count the phenomenology of the H;Sbc emitting region, as well 
as inferences about BLR structure are quite similar over a 5 dex 
luminosity range that includes some of the most luminous known 



Marziani et al.: H/3 in Intermediate z Quasars: III. Protile Analysis 



11 



quasars. This overall scenario is described in more detail by the 
following results of the present paper: 



1. A minimum requirement for quantifying and interpreting 
BLR properties in Type 1 AGN involves the Pop. A-B di- 
chotomy. Spectroscopic averaging without such a discrimi- 
nation appears to be equivalent to discussing stellar proper- 
ties without consideration of the OBAFGKM spectroscopic 
sequence 

2. Pop. A sources show a minimum FWHM Hp that increases 
with source luminosity from ~500-1000 km s"' to 3000- 
3500 km s The best boundary to this lower envelope fa- 
vors an exponent a ^ 0.67 in the Kaspi relation. Pop. A 
sources span a FWHM H/J range of ~ 4000 km s ' driven 
by source orientation and virial motions in an accretion 
disk. The virial assumption implies a Mbh range of 4dex 
(log Mbh= 6-10) for this high accreting population. In the 
ISAAC sample studied in this paper Pop. A sources show 
-0.7 < log Lboi/iEdd^ 0.18. 

3. HjS emission in Pop. B sources involves both unshifted broad 
(HjSbc) and redshifted very broad line (HjSvbc) components. 
Pop. B shows a "super" Baldwin effect with W(Hy6Bc) de- 
creasing by a factor 4 with increasing Lboi- At the same time 
W(I-l;6vBc) remains almost constant implying that it arises 
in a photo-ionized optically thick medium. FWHM(Hj6bc) 
component shows no luminosity dependence after correction 
for HjSvBC broadening. This yields an Mbh range similar to 
that for Pop. A. Pop. B is a lower accreting population with 
-1.1 $ log Lboi/i-Edd^ -0.2 for the objects of this paper 
One should bear in mind the caveats of ^4.61 in analyzing 
the Lboi/^Edd hmits, and that 80% of Pop. B sources of the 
ISAAC sample show log Lboi/i'Edd~ -0.4. Considering the 
blurring by errors in Lboi/^'Edd estimates, and results of pre- 
vious papers, intrinsic Lboi/i'Edd ranges could be k 0.01 - 
0.2 (Pop. B) and a; 0.2 - 1 (Pop. A), with a separation at 
^boi/^Edd~ 0.2 + 0.1. Pop. B sources in the SDSS sample 
seem to remain a population of lower Eddington ratio even if 
the extreme correction to FWHM implied by removal of the 
HjSvBC is applied ( d1 ). 

4. Most Type 1 AGN do not show broad line profiles consistent 
with simple accretion disk models. Perhaps Lorentzian Pop. 
A profiles can be accommodated with thin/slim disk emis- 
sion models (if disk emission is very extended) but Pop. B 
cannot. The redshifted Hy6vBC in Pop. B sources can not be 
considered part of a double peak signature since the asym- 
metry and the centroid shift at j intensity (near the pro- 
file base is usually too large to be attributed to a gravita- 
tional redshift. Arp 102B - the prototypical "double-peaked" 
Babner line emitter - shows a redshift at Une base that is 
modest, and consistent with the value expected for gravita- 
tional redshift at the inner emit ting radius o f the accretion 
disk according to the model of IChen et al.l (Il989l) . This is 
however not the case for many Pop. B objects that have 
been included in so-called double-peaked s ource compila- 
tions (i.e. as accretion disk candidates: e.g., IStrateva et al.L 
I2003h . 

5. HjSvBC is the dominant broad line component at high lumi- 
nosity (from ~ 25% to 60% of the H/? flux over our Lboi 
range) for Pop. B sources. The geometry and kinematics of 
the VBLR are yet unclear The HySvBC redshift Av ^ 1-2000 
km s"' is unlikely to be gravitational in nature. The increase 
in dominance with redshift favors a connection to the hy- 
pothesized infall involved with the Mbh growth. Another 



possibility inv olves photon downshifting via some form of 
scattering (e.g. lLaoAl2006l and references therein). 
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Fig. 1. Calibrated VLT-ISAAC spectra for 30 new intermediate-redshift quasars. Abscissae are rest-frame wavelength in A, ordinates 
are rest-frame specific flux in units of 10 erg s ' cm ' A The spectrum behavior of HE 0035-2853, HE0043-2300, HE 
0058-3231 is fully unreliable beyond ^ 5250 A due to heavy atmospheric absorption also affecting the calibration standard. 
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Fig. 2. Continuum-subtracted spectra. Left-hand panels: the H/3 spectral regions. The thin green lines show the Feii emission. Right- 
hand panels: spectra after continuum and Feii subtraction. The HyS profile (broad + very broad) is marked by the thick line, blue on 
the short-wavelength side with respect to the rest wavelength, red on the long-wavelength side. Abscissa is rest frame wavelength 
in A, ordinate is specific flux in units of 10 erg s ' cm"^ A 
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koi [er'gs s-i] [ergs s"'] 

Fig. 3. FWHM(Hy6) in km s"' (a). Asymmetry Index (b), centroid at 3/4 fractional heights, in km s"' (c), and centroid at 1/4 
fractional height, in km s~' (d ) as a function of bolometric luminosity in units of erg s"'. Small gray circles: data from the SDSS 
sample of IZamfir et alj (|2008|) : larger filled circles: ISAAC sample. Circled data points identify RL sources. In the upper left panel, 
the continuous line shows an unweighted Isq best fit. The dotted lines identify the +2cr confidence limit for A.I., and centroids equal 
to zero. 
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Fig. 4. The optical plane of Eigenvecto r 1 with the VLT sources presented in this paper and in Paper 1 and II (large red filled circles), 
and the SDSS sample of IZamfir et al.l (l2008i small grey circles)). Circled data points identify RL sour ces. Abscissa is Fen/ 14570 
prominence /Jpen, ordinate is FWHM(Hj8) in km s ' . The plane is subdivided in spectral types following lSulentic et al.l (|2002|) . 
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Fig.5. Median spectra for bins Al, Bl, A2, B2, as defined in Fig.|4] Original spectral have been normalized to unity si A - 5100 A. 
The adopted continuum is shown as a solid line. 
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Fig. 6. Analysis of the HyS profiles for spectral types Al, A2, Bl, B2. Median spectra are shown after continuum subtraction (thin 
line); the model fit obtained with specfit is shown as a thick dashed magenta line. The H/3 profile is represented either by a 
Lorentzian function or by the sum of two Gaussians (not shown individually; thick lines). The thin green line shows the Feiiopt 
contributions, and the yellow lines the [Oiii]/l/l4959,5007 and H/3nc emission. Residuals of the fitting procedure are displayed in 
the window immediately below the one showing the spectral components. Abscissa scale is rest frame wavelength in A (top), and 
radial velocity in km s"' from rest wavelength of Hfi (bottom). 
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Fig. 7. Luminosity effects in the median spectra of Pop. A and Pop. B sources, binned in decades of luminosity in the luminosity 
ranges 43 < logLbo/ < 49. Abscissa is rest frame wavelength in A, ordinate is normalized intensity (set to unity at /I = 5100 A). 
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Fig. 8. Median spectra analysis of the Hj8 profiles as a function of luminosity, for Pop. A and B (right panels). The median spectra 
have been continuum-subtracted. Lines follow the same coding of Fig. |6] Residuals of the fitting procedure are displayed in the 
window immediately below the one showing the spectral components. 
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Fig. 8. Cent. 
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Fig. 9. Distribution of unshifted and shifted H/J line centroids at 1/4 fractional intensity as a function of the bolometric luminosity 
(logLboh binned over AlogLboi - 1, with the naming convention of Tables |6] |7] [SJ for the Pop. B sources of the merged SDSS 
ISAAC sample. Dark grey (bottom): unshifted; grey (middle): redshifted; pale grey (top): blueshifted. "Unshifted" means that the 
centroid measure is within a 2cr confidence range of +400 km s"' . 
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Fig. 10. Left: distribution of the ISAAC sources whose H/?bc can be fit by a Lorentzian (dashed histogram) and by a Gaussian as a 
function of bolometric luminosity in ergs s Right: the same ISAAC sources in the FWHM(HyS) vs. bolometric luminosity plane. 
Each source is identified according to the best-fitting function for the core of Wfi: Lorentzian (L), Gaussian (G), or uncertain (U, 
undecided). The lower dot-dashed line sets the limit of Pop. A at 4000 km s the upper line traces the luminosity-dependent limit 
that may be more appropriate. In between the two lines, an intermediate population of sources (called Pop. M) show predominantly 
Lorentzian profiles. 



26 



Marziani et al.: in Intermediate z Quasars: III. Protile Analysis 



o 
o 
o 

CD 



O 
O 
O 
C\2 



o 



^ 2 




o — 



42 



Fig. 11 . FWHM(H8) in km s 
(|2008|) with the "deep limit 



44 



46 



48 



' as a function of bolometric luminosity. Small gray circles: data from the SDSS sample oflZamfi r et al 
< 17.5; larger filled circles: ISA AC sample; blue filled tr iangles: data for NLSyls from iZhou et al 



( 2006 ). The open circles are the data points from the ATLAS of iMarziani et alj (l2003al) . The horiz ontal straight, dot dash ed lines 
show the limit of classical NLSyls sources (2000 km s ') and the limit derived for Population A of ISulentic et al.l (l2000ah at low z 
i.e., 4000 km s ' . The continuous line show the minimum FWHM expe cted in case of viria l motion, for tw o values of the pow er-law 
index in the relationship between BLR size and luminosity: a = 0.52 dBentz et al.L |2006|) and a = 0.67 ("K aspi et"an , |2005', Paper 
II). The thick dot-dashed curve at FWHM ^ 4000 km s ' traces the luminosity-depedent boundary between Pop. A and B partly 
adopted in this paper. The red thick lines show the FWHM increase in luminosity expected for Lboi/i'Edd= 0.15 (assuming a = 0.67). 
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Table 1. Basic Properties of Sources and Log of Observations 
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1.10-1.04 


10 


HE0940-1050 




16.96 


3.0932 




I 


-29.82 


< 0.10 


2005/06/10 


sK 


180 


24 


1.18-1.50 


10-45 


HE 1 039-0724 




17.16 


1.4584 


0.0002 


1,2,3 


-27.86 


< 0.20 


2005/06/23 


J 


180 


20 


1.19-1.44 


20 


HE 1 120+0154 


UM 425 


16.31 


1.4720 


0.0004 


1,2,3 


-28.70 


-0.57 


2005/06/24 


J 


180 


16 


1.42-1.75 


35 


HEI347-2457 




16.83 


2.5986 




I 


-29.64 


-0.68 


2005/04/15 


sH 


180 


28 


1.16-1.54 


65-10 


T TT^ -t ^ A a , f\ t to 

HEI348+01 18 


Q 1348+0118 


17.08 


1.0873 


0.0012 


I 


-27.05 


< 0.24 


2005/06/13 


sZ 


180 


12 


1.20-1.29 


15 


HEI349+0007 


UM617 


16.83 


1.4442 


0.0010 


1,2,3 


-28.00 


-0.18 


2005/06/24 


J 


180 


20 


1.25-1.52 


15 


HEI409+0101 




16.92 


1.6497 


0.0039 


I 


-28.30 


0.40 


2005/07/21 


J 


180 


20 


1.18-1.38 


55-25 


HEI430-004I 


Q I430-004I 


16.72 


I.I2I6 




I 


-27.60 


-0.39 


2005/07/21 


sZ 


180 


16 


1.43-1.78 


20 


HE 1 505 +02 1 2 




17.06 


1.0943 




I 


-27.28 


0.51 


2004/07/18 


sZ 


180 


16 


1.26-1.53 


50 


HE 1 505 +02 1 2 






1.0930 


0.0240 


I 






2005/05/21 


sZ 


180 


20 


1.28-1.57 


30 


HE2147-3212 




16.84 


1.5432 




I 


-28.15 


< 0.14 


2004/07/05 


J 


180 


20 


1.06-1.19 


15 


HE2 156-4020 




17.39 


2.5431 


0.0010 


I 


-28.75 


-0.09 


2004/07/10 


sH 


180 


32 


1.05-1.20 


15-5 


HE2335-3029 




16.96 


1. 1194 


0.0016 


I 


-27.21 


<0.20 


2004/07/1 1 


sZ 


180 


16 


1.58-1.31 


25 



Apparent B lohnson magnitude corrected because of Galactic absorption. The value of HE 0436-3709 referes to the b Stromgren magnitude. 
^ Redshift, with uncertainty in parentheses. 
Lines used for redshift calculations: 1: 11^6, 2: [Oiii]/15007. 3: [Oiii]14959. 

Absolute B magnitude, computed for //o = 70 km s"' Mpc"', Qm = 0.3, f^A = 0.7, and fc-correction spectral index a = 0.6. 
° Decimal logarithm of the specific flux ratio at 6cm and 4400 A (effective wavelength of the B band). Upper limits are from the NVSS (~ 2.5 

mjy), and would place all undetected sources in the RQ domain. 
^ Date refers to time at start of exposure. 
^ Photometric band. 

^ Detector Integration Time (DIT) of ISAAC, in seconds. 
' Number of exposures with integration time equal to DIT. 
j Airmass at start and end of exposure. 

S /N al continuum level in the proximity of H/3. Two values are reported in case of different S/N on the blue and red side of H/S (blue side 
first). The 5 /N value is with N estimated at a 2cr confidence level i.e., 2 times the rms. 
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Table 2. Measurements of Fluxes, Equivalent Widths and FWHM of Strongest Lines 



Object name 




W(H/3)'' 


F(Feii/l4570)'^ 


FWRMCFen)"* 


(1) 


(2) 


(3) 


(4) 


(5) 



HE0035-2853 
HE0043-2300 
HE0058-323I 
HE0I05-3006 
HE0I09-3518 
HE0I19-3325 
HE0142-4619 
HE0203-4627 
HE0205-3756 
HE0224-3826 
HE0229-5603 
HE0239-4809 
HE0248-3628 
HE025I-5550 
HE0349-5249 
HE0359-3959 
HE0435-4312 
HE0436-3709 
HE0940-1050 
HE1039-0724 
HE II 20+0 1 54 
HE 1347-2457 
HE1348+0I18 
HE 1349+0007 
HE1409+0I0I 
HE1430-004I 
HE1505+0212 
HE2I47-3212 
HE2I56-4020 
HE2335-3029 



115 +3° 
214 +i 
173 

58 
831 

93 
274 
159 
343 
123 
118 
111 
250 
353 
289 



60+" 



+28 
+Foo 

+if 

-30 
+20 

d8 

-40 
+ 15 

rf8 



+40 
-40 
+50 

+-r? 



165 
139 

Q99 +200 
-100 

138 
270 

264 +1 
611 +*'^ 

-60 

129 

44^ +r50 
^^-^ -100 

87 i 

176!^° 

+20 



151 
212 
142 



;?8 



74+10 

100 +1^ 

83 !|f 

119!!^ 

133 

77 +'2 

74^1? 
-11 

77 

' ' -23 
77 +10 

44+ff 

40+"' 

71!!1 
63!" 

45 t 
118!?^ 

43 f 
75 
63 



64 

95 +" 

^-^ -13 

76 

-14 

91 

_22 

50 1 
59 

82+1^ 



+ 10 

.-.■to 



-40 
hlO 

+'Po 

-20 
10 i 



82 
72 



lOI _ 

82!fo 
60 !| 

39 !f 

95 

64 

104 

76 

102 

98 

5"° _2n 
04 +10 
-10 

102 

119+1§ 

67 

186 2° 
III 
86 

505 !!® = 
26 +f, 
35 

-15 

252 
255 ^^^ 
47 

^' -10 
979 +20 

7Q +30 

-10 

350 

-150 

139 

-40 

95 !'^ 
25 



6400 
3900 
5800 
6000 

4500 
4100 
6600 
4900 
5100 
3000 
2500 
4300 
2200 
4700 
2600 
5900 
5200 
3000 



+2000 

m 

-1400 
+2000 

m 

-1800 
700 



+ 12O1 
- 1 ''Oi 



-600 
+600 
-2200 
+800 

;l588 



-900 



5100 +™ 
4300 

t JVJU 2400 

4000 +^§00 

6300 

1000 +J™ 

2700!™ 

5800!™ 

5200 

2000 +™ 



Rest frame flux of IIy6 in A ±2cr confidence level uncertainty. 
Rest frame equivalent width of HyS in A ±2cr confidence level uncertainty. 
Rest frame flux of the Feii/i4570 blend in units of erg s"' cm"^. 
FWHM of lines in the blend. See text for details. 

Rest frame flux of the Feii/i4570 blend estimated from the Fen/i5270 blend. See text for details. 
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Table 3. Measurements of Fluxes and Equivalent Widths of Narrow Lines 



Object name 




WdiSNc)" 


F([Oin]^5007)'^ 


W([Oiii]/i5007)'' 


(1) 


(2) 


(3) 


(4) 


(5) 



12 4 

i^-t -2.0 

27.8 

13 1 
1 2 

-'■^ -1.2 

247 5 +25.0 
1 6 

-1.0 

20 

8-3!:^ 

33.0!]^:° 
<3.0 
4 

-4.0 
7 9 +2.0 
'•^ -2.0 

15 +^-' 
49 

0.8-^ 
13 4 " 

T 9 +1.0 
•'■^ -1.0 
1 I 1 +100.0 
^^■^ -11.0 

26 2 
52 6 +'* ° 
<3.0! 
60 7 

-20.0 
Q 4 +2.0 
^■^ -2.0 

35 7 

-10.0 

<2.0 
<2.0 

i.8![:« 

3.8 +Jj 

19.0 +tn 



HE0035-2853 
HE0043-2300 
HE0058-3231 
HE0105-3006 
HEO 109-35 18 
HEOl 19-3325 
HE0142-4619 
HE0203-4627 
HE0205-3756 
HE0224-3826 
HE0229-5603 
HE0239-4809 
HE0248-3628 
HE025 1-5550 
HE0349-5249 
HE0359-3959 
HE0435-4312 
HE0436-3709 
HE0940-1050 
HE1039-0724 
HEl 120+0154 
HE 1347-2457 
HE1348+0118 
HE 1349+0007 
HE1409+0101 
HE1430-0041 
HE1505+0212 
HE2147-3212 
HE2156-4020 
HE2335-3029 



4 0+'" 

-1.0 

12.0 ^Ifo 

0-3 if 

74 4 +10.0 

0-8 1? 
T q +7.0 
-'■'^ -2.0 

1 +-" 

-1.0 

<7.0 
1 0+'" 

-1.0 

< 1.0 : 
2-8!?:^ 

4 8 +-■" 

-5.0 

3 3 +^-'' 
4^1:8 
0.4 J 

1 5 

5:0 f 

9 9 +i:§ 

5-0 

18.7 

2 6 +"■' 
34 8 +^-'' 

-15.0 

<2.0 

< 1.0 

2 0+'" 
-1.0 

5.0 !2 " 

4 

-2.0 



-4.2 
+1.7 



-0.4 
1.5 



2.1 
4.1 
1.0 
0.4 

3.5 _ 

7 +'-^ 

^•^ -1.0 

5 +'■" 
"•-^ -0.5 

< 1.5 
+° '' 

<0.7 
1.1 
0.8 
1.1 
0.7 
2.3 
0.3 
0.9 
0.3 
2.7 
0.5 
1.2 
2.9 
1.5 
7.1 



+0.2 
-0.4 
+0.3 



-0.1 
+0.3 

-1.1 
+0.3 

+8:9 
-3.2 
< 1.2 

<0.3 

1.1 
1.9 



2.0 



6.7 !| ? 
10.1 
7.9 *\\ 

1.8 i^ 
38.5 +" 

1- 6 -If 

4.2 +" 

7 4 +^-^ 

< 1.9 

2- 7 1-^ 
9 Q +o;8 
^■^ -0.8 
2 6 +" * 

-OA 

10.3 

2 5 +"-^ 
-1.1 

"•^ -oj 
10.1 +^0 

2 1 +"■' 
^•i -0,7 

6 ' 

15 5 +^-' 

iJ-J -4.1 

13.1 +2-3 

c Q +0.6 
^■^ -0,6 
10.1+]^ 

6.1 !|i 
7.8 

< 1.2 
<0.7 
1.1 +"■* 

1 5 +^-^ 

-'A 

10-2!?;: 



" Rest frame flux of H/Jnc in units of 10 '^ erg s ' cm ^. 
^ Rest frame equivalent width of H/Jnc in A. 

Rest frame flux of [Oni]/15007 in units of 10"" erg s"' cm"^. 
'' Rest frame equivalent width of the [Om]/i5007 line in A. 
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Table 4. Hfi Line Profile Measurements 



Source 


FWZP 


■ h 


FJ J\ Til 


■ h 


A.I.*^ 


Ab 
A 


Kurt. 


A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HE0035-2853 


23300 


3600 


7400 


500 


0.19 


0.07 


0.37 


0.04 


HE0043-2300 


25000 


4500 


4600 


300 


-0.03 


0.10 


0.33 


0.06 


HE005 8-3231 


25200 


4400 


5800 


500 


0.22 


0.08 


0.30 


0.04 


HEO 1 05-3006 


20700 


5000 


5200 


400 


0.10 


0.17 


0.28 


0.08 


HE0I09-3518 


24400 


5900 


4000 


300 


-0.01 


0.10 


0.33 


0.06 


HEO 11 9-3325 


25300 


6800 


4500 


300 


0.03 


0.09 


0.41 


0.06 


HE0142-4619 


26900 


4200 


3600 


300 


0.24 


0.17 


0.26 


0.08 


HE0203-4627 


27800 


5100 


7200 


500 


0.12 


0.09 


0.32 


0.05 


HE0205-3756 


20700 


4400 


4400 


400 


0.41 


0.06 


0.16 


0.03 


HE0224-3826 


20800 


2900 


6500 


400 


0.10 


0.07 


0.40 


0.05 


HE0229-5603 


26500 


6300 


4000 


300 


0.03 


0.10 


0.31 


0.06 


HE0239-4809 


19700 


3900 


4700 


400 


0.19 


0.14 


0.27 


0.06 


HE0248-3628 


22200 


4400 


3900 


300 


-0.19 


0.10 


0.30 


0.05 


HE025 1-5550 


25700 


4400 


4500 


400 


0.47 


0.07 


0.19 


0.03 


HE0349-5249 


26500 


4000 


5900 


800 


0.37 


0.07 


0.27 


0.04 


HE0359-3959 


18500 


4800 


4300 


400 


-0.15 


0.09 


0.30 


0.05 


HE0435-4312 


25300 


6100 


5500 


300 


0.08 


0.14 


0.38 


0.08 


HE0436-3709 


26700 


4300 


6600 


500 


0.26 


0.09 


0.32 


0.05 


HE0940-1050 


15400 


1400 


3500 


200 


0.53 


0.04 


0.22 


0.02 


HE1039-0724 


22900 


3300 


11300 


500 


0.63 


0.04 


0.48 


0.03 


HEl 120+0154 


25800 


4300 


6900 


500 


0.33 


0.07 


0.35 


0.04 


HE 1347-2457 


24700 


3600 


6800 


500 


0.20 


0.09 


0.31 


0.05 




ZjOUU 




J /UU 


ZUU 


U.Uo 


n 1 


n 

yj.JO 


U.Uo 


HE 1349+0007 


23400 


3500 


6200 


400 


0.27 


0.10 


0.34 


0.06 


HE1409+0101 


21600 


4600 


8300 


900 


0.48 


0.07 


0.28 


0.03 


HE 1430-0041 


19100 


3600 


4700 


300 


-0.05 


0.10 


0.33 


0.06 


HE1505+0212 


18400 


2300 


5200 


200 


-0.02 


0.25 


0.37 


0.11 


HE2147-3212 


25900 


3900 


5800 


400 


0.01 


0.10 


0.33 


0.05 


HE2156-4020 


26200 


4400 


5000 


500 


0.27 


0.10 


0.26 


0.05 


HE2335-3029 


18100 


4000 


3100 


200 


0.04 


0.11 


0.38 


0.07 



" In units of km . 

^ 2cr confidence level uncertainty. 

Asymmetry index defined as in Marziani et al. (1996). 

Kurtosis parameter as in Marziani et al. (1996). 
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Table 5. Hfi Line Centroids at Different Fractional Heights 



Source 


c(0/4) 


A li b 


C(i/4j 


A a b 
A 


c(l/2) 


Aab 


c(3/4) 


Aab 
A 


C(U.9) 


Aab 

A • 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


HE0035-2853 


3100 


1800 


1900 


400 


930 


270 


490 


200 


370 


130 


HE0043-2300 


100 


2200 


100 


400 


130 


160 


130 


120 


130 


80 


HE0058-3231 


3500 


2200 


2500 


500 


1010 


240 


710 


150 


640 


100 


HE0105-3006 


2900 


2500 


1000 


1000 


290 


180 


250 


140 


240 


90 


HE0109-3518 


400 


3000 





300 


40 


140 


40 


110 


40 


70 


HE0119-3325 


800 


3400 


100 


300 


40 


130 


30 


120 


20 


80 


HEO 142-46 19 


1300 


3400 


100 


300 


40 


130 


30 


120 


20 


80 


HE0203-4627 


-300 


2100 


1000 


800 


200 


150 


110 


100 


90 


60 


HE0205-3756 


2900 


2200 


2600 


500 


580 


220 


80 


110 


10 


60 


HE0224-3826 


3400 


1500 


1100 


400 


650 


200 


400 


180 


300 


120 


HE0229-5603 


1500 


3200 


-100 


300 


-110 


150 


-110 


110 


-110 


60 


HE0239-4809 


-100 


2000 


200 


600 


-330 


180 


-330 


130 


-310 


70 


HE0248-3628 


500 


2200 


-400 


400 


10 


160 


130 


110 


150 


60 


HE025 1-5550 


2900 


2200 


2900 


600 


90 


200 


-60 


120 


-90 


70 


HE0349-5249 


2700 


2000 


2200 


500 


640 


420 


100 


150 


20 


90 


HE0359-3959 


-700 


2400 


-800 


400 


-320 


180 


-140 


120 


-90 


70 


HE0435-4312 


2700 


3100 


400 


600 


100 


160 


50 


150 


40 


100 


HE0436-3709 


2600 


2100 


2100 


600 


570 


240 


310 


180 


240 


120 


HE0940-1050 


3200 


700 


2700 


200 





90 





90 





70 


HE 1039-0724 


2000 


1600 


1900 


300 


1270 


250 


140 


830 


-1380 


160 


HE1120-H0154 


2700 


2200 


1300 


500 


60 


260 


-260 


190 


-350 


120 


HE1347-2457 


2400 


1800 





600 


-390 


270 


-570 


190 


-610 


110 


rlj_/l jH-oi^U i lo 




ZZUU 


^00 


zLOO 


70 


1 1 n 

1 iU 


^0 


1 00 
lUU 




/U 


HE1349+0007 


3000 


1700 


1400 


700 


180 


200 


-10 


170 


-60 


110 


HE1409+0101 


800 


2300 


1600 


500 


510 


460 


-700 


210 


-860 


110 


HE1430-0041 


-400 


1800 


200 


400 


200 


170 


200 


130 


200 


80 


HE 1505 +02 12 


-800 


1200 





900 


-390 


110 


-60 


260 


30 


80 


HE2147-3212 


-1300 


2000 


-100 


500 


-70 


210 


-70 


160 


-70 


100 


HE2 156-4020 


1700 


2200 


1600 


500 


390 


250 


150 


130 


100 


80 


HE2335-3029 


800 


2000 


100 


300 


10 


90 


-10 


80 


-10 


60 



^ In units of km s ' . 

^ 2cr confidence level uncertainty. 
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Table 6. Measurements on the Broad Lines of Median Spectra 



Object name 


(2) 


W(Feia4570)" 


FWHM(Feii)'' 


F(H8bc)/F(H8)'-- 


FWHM(Hy3Bc)'' 


log Mbh Wnc 
(7) 


logLboi/i'Edd H/3bc^ 


Al 


72+11 


26 +3 


2700+!!°° 


1.00 








A2 
Bl 


65 

86!" 


49 +" 

-11 

26 +j 


3700 +^o°o 
5200 +^'*"° 


1.00 
0.27 


4000 






B2 


44 +^ 

-14 


5000 +^" 


0.32 


4000 






A 
M 


61 

67+'" 


25 +^ 
35 !' 


2700 !|5°SJ 
3800 


1.00 
1.00 








MB 


86 


31 !° 


5000 !jf 


0.27 


4100 






43A 
44A 
45A 
46A 
47A 


91 

-20 

69 

86 *\l 

68!!? 
60 _„ 


36 +' 
38!1« 

+10 
^■^ -10 

47!|« 
301 

27 _g 


3000 +„!; 
2600 !5"» 
2800 !J58 
3000 I™ 
3000 

Qonn +1500 
3800 _i, 00 


1.00 
1.00 
1.00 
1.00 
1.00 

1 nn 


... 


6.1 
6.8 
7.8 
8.6 
9.6 

1 n Q 


-0.74 

-0.47 

-0.43 

-0.26 

-0.20 
I n 1 1 


43B 


130!^° 


8!j« 


f 


0.59 


4600 


7.1 


-0.68 


44B 
45B 


125 

111 

ill 20 


38 !2o 
29 +5 

-15 


5600 +''™ 
4900 +''"' 


0.49 
0.35 


4700 
4400 


7.7 
8.4 


-1.37 
-0.98 


46B 
47B 
48B 


QT +10 

92+13 

-14 

75 +' 


22 +' 

-10 

38 +^ 
12 1 


5900 +'^'' 

J^i-w J 200 

4900 +"''"' 

t7UU 2000 

4600 +'-'"' 


0.37 
0.27 
0.23 


4800 
4000 
4300 


9.1 
9.6 
10.3 


-0.73 
-0.24 
+0.03 



" Equivalent width of US (Hj8bc+ HjSvbc) and Feii/I4570in A ±2cr confidence level uncertainty. Note that those values have been computed on 
median spectra with flux normalized to unity al A = 5100 A. Considering that the continuum shape is not flat, but that there is however little 
dispersion in continuum shape across the median spectra, it is W(H/3) a I(II/?)/l.l, and W(Feii/i4570)a;I(Feii/i4570)/1.25 A. 

^ FWHM of lines in the blend in units of km s"' computed by specf it as for the individual sources. Uncertainty is at ±2cr confidence level. 
See text for details. 

Intensity ratio of the H/Jbc to total H/J line emission i.e., H/Jbc and H/Jvbc- 
FWHM of the H/3bc component i.e., after removing H/3vbc. 

Logarithm of Mbh, in solar masses, and of Lboi/^^Edd- Values have been computed following Paper II, using the FWIIM(Hj8bc) reported in 
Col. (6), and assuming the average bin luminsosity. Values are therefore only indicative. No Mbh or Lboi/^Edd has been computed for median 
in spectral types since they are normalized median spectra made regardless of their luminosity. 

Feiiop, too faint for FWHM to be meaningfully constrained. 
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Table 7. Hj6 Line Profile Measurements on Spectral Types 



Source 


rWZl 


A a b 

A 




A ii b 

A 


A.l. 


A 


Kurt."* 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A 1 




IH-UUU 






U.U i 


U. iU 


fl 

U.J J 


fl flfi 


A2 


41100 


16000 


4100 


300 


0.05 


0.10 


0.33 


0.06 


Bl 


34700 


6000 


5900 


600 


0.37 


0.10 


0.25 


0.04 


B2 


29200 


5000 


5700 


500 


0.26 


0.09 


0.29 


0.05 


A 


37100 


15000 


3700 


300 


0.02 


0.10 


0.33 


0.06 


M 


43700 


17000 


4400 


300 


0.01 


0.10 


0.33 


0.06 


BM 


34400 


6000 


6300 


700 




fl flQ 




fl fl4 


43A 


19800 


8000 


2000 


100 


-0.05 


0.11 


0.33 


0.07 


44A 


20700 


8000 


2100 


100 


-0.07 


0.11 


0.33 


0.07 


45A 


27300 


11000 


2700 


200 


-0.01 


0.10 


0.33 


0.06 


46A 


31000 


12000 


3100 


200 


-0.04 


0.10 


0.33 


0.06 


47A 


40000 


16000 


4000 


300 


0.01 


0.10 


0.33 


0.06 


48A 


40900 


16000 


4100 


300 


0.05 


0.10 


0.33 


0.06 


43B 


27900 


6400 


5200 


300 


0.06 


0.09 


0.41 


0.05 


44B 


31600 


6600 


5700 


300 


0.07 


0.11 


0.39 


0.06 


45B 


34800 


6400 


6000 


400 


0.20 


0.13 


0.30 


0.06 


46B 


32900 


5900 


6500 


500 


0.22 


0.10 


0.31 


0.05 


47B 


36900 


6400 


5900 


500 


0.29 


0.11 


0.25 


0.05 


48B 


29400 


4600 


9300 


900 


0.44 


0.06 


0.29 


0.03 



" In units of km s"' . 

2cr confidence level uncertainty. 

Asymmetry index defined as in Marziani et al. (1996). 

Kurtosis parameter as in Marziani et al. (1996). 



Table 8. H^S Line Centroids at Different Fractional Heights on Spectral Types 



Source 


c(0/4)" 


^a.b 


c(l/4)" 


^a,b 


c(l/2)^ 


^a.b 


c(3/4)" 


^a.b 


c(0.9)" 


^a.b 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


Al 





6900 





300 


-30 


130 


-30 


100 


-30 


60 


A2 


-100 


8200 


-200 


300 


-170 


150 


-170 


110 


-170 


70 


Bl 


2400 


3000 


2100 


700 


130 


280 


-130 


150 


-170 


100 


B2 


2200 


2600 


1700 


600 


380 


230 


170 


150 


130 


100 


A 


-100 


7400 


-100 


300 


-70 


130 


-70 


100 


-70 


60 


M 





8700 





400 


-40 


160 


-40 


120 


-40 


70 


BM 


2400 


2900 


2200 


700 


270 


330 


-60 


160 


-120 


100 


43A 


100 


3900 


100 


200 


90 


70 


90 


50 


90 


30 


44A 


200 


4100 


100 


200 


130 


70 


130 


60 


130 


40 


45A 





5400 





200 


30 


100 


30 


80 


30 


50 


46A 


100 


6200 


100 


300 


120 


110 


120 


90 


110 


50 


47A 





7900 





300 


-20 


140 


-20 


110 


-20 


70 


48A 


-100 


8100 


-200 


300 


-170 


150 


-170 


110 


-170 


70 


43B 


2100 


3200 


200 


300 


40 


150 





140 


-10 


100 


44B 


1200 


3300 


200 


500 


10 


170 


-30 


150 


-40 


100 


45B 


1700 


3200 


1200 


800 


140 


210 


40 


160 


20 


100 


46B 


2400 


3000 


1600 


700 


370 


230 


190 


170 


150 


110 


47B 


2100 


3200 


2000 


800 


220 


260 


30 


160 





100 


48B 


3100 


2300 


2700 


400 


1550 


440 


-20 


240 


-230 


120 



" In units of km s ' . 

2cr confidence level uncertainty. 




Fig. 12. Median spectra for Population A, intermediate spectra in region M, and modified Population B (BM), as defined in Fig.fTOl 
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Fig. 13. Median spectra analysis of the H/3 profiles for Population A, "middle" Population M, and Population MB. Units and line 
formats are the same as for Fig. [8] 
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Fig. 14. Distribution of sources whose HySec is of Lorentzian (hatched histogram) or Gaussian type as a function of black hole mass 
in units of solar masses (left) and of Eddington ratio (right). 
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Fig. 15. Distribution of the ISAAC sources in the Lboi/i^Edd vs. Mbh (in solar units) plane. Each source is labeled as in Fig.fTOl 



Marziani et al.: H/3 in Intermediate z Quasars: III. Profile Analysis, Online Material p 12 



4000 



3000 - 



-SOOO 



6 

^000 



■1000 



m 



9 



_L 



_L 



-0.5 



0.5 



4000 



3000 




1000 



9 9.5 10 



10.5 



0.8 p- 
0.7 r 
0.6 - 
0.5 r 
0.4 ~- 
0.3 ^ 
, 0.2 - 
0.1 
\^ 
-0.1 ^ 
-0.2 ~- 
-0.3 ^ 
-0.4 - 
-0.5 



m 



® 



i: 



_L 



_L 



■1 



-0.5 
logio (L/LEdd) 



0.5 




9.5 



10.5 



Fig. 16. Dependence on black hole Mbh and Eddington ratio for HyS c(|) (upper panels) and A.I. (lower panels). c(|) is in units 
of km s Mbh in units of solar masses. Circled data points identify RL sources. Typical error bars at a 2cr confindence levels are 
displayed in the corners of each panel. Significant correlations are presents for c{j) and A.I. vs. log Mbh; the thick lines are best fits 
with the unweighted least squares method. The dashed lines represent unweighted least squares fits restricted to Pop. B. 
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Fig. 17. Hy8 median profiles for sources in 2 ranges of Mbh and of Lboi/iEdd- The upper panels show the higher Lboi/iEddfange 
(log Lboi/i'Edd^ -0.5), while the bottom panels are for -1.1 ^ log Lboi/^Edd< -0.5. The lower mass range in the left panel is 
8.5 ^ logMBH< 9.5. Meaning of lines is the same of Fig. |6]and[8] See text for more details. 



